The Effect of a Chameleon Scalar Field on the Cosmic Microwave Background 
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We show that a direct coupling between a chameleon-like scalar field and photons can give rise to a 
modified Sunyaev-Zel'dovich (SZ) effect in the Cosmic Microwave Background (CMB). The coupling 
induces a mixing between chameleon particles and the CMB photons when they pass through the 
magnetic field of a galaxy cluster. Both the intensity and the polarization of the radiation are 
modified. The degree of modification depends strongly on the properties of the galaxy cluster such 
as magnetic field strength and electron number density. Existing SZ measurements of the Coma 
cluster enable us to place constraints on the photon-chameleon coupling. The constrained conversion 
probability in the cluster is 7 , coma(204 GHz) < 6.2 x 10~ 5 at 95% confidence, corresponding to an 



upper bound on the coupling strength of t^j ell) < 2.2 x 10 -8 GeV -1 or g^ olmo> < (7.2 - 32.5) x 



(Kolmo) 



10" 



' GeV , depending on the model that is assumed for the cluster magnetic field structure. We 



predict the radial profile of the chameleonic CMB intensity decrement. We find that the chameleon 
effect extends further towards the edges of the cluster than the thermal SZ effect. Thus we might see 
a discrepancy between the X-ray emission data and the observed SZ intensity decrement. We further 
predict the expected change to the CMB polarization arising from the existence of a chameleon-like 
scalar field. These predictions could be verified or constrained by future CMB experiments. 



PACS numbers: 



I. INTRODUCTION 



Measurements of the current cosmic acceleration of the 
universe have led to theories of an unknown energy com- 
ponent with negative pressure in the Universe. This dark 
energy is generally modelled by a scalar field rolling down 
a flat potential. On Solar-system scales this scalar field 
should appear effectively massless. The origin of the dark 
energy scalar field is still highly speculative; some sug- 
gestions have involved the moduli fields of string theory. 
However, if a near-massless scalar field existed on Earth 
it should have been detected as a fifth force unless the 
coupling to normal matter was unnaturally suppressed. 
Khoury and Weltman [l[ introduced the chameleon scalar 
field which has gravitational strength or greater coupling 
to normal matter but a mass which depends on the den- 
sity of surrounding matter. On cosmological scales the 
chameleon successfully acts as a low-mass dark energy 
candidate, while in the laboratory its large mass allows 
it to evade detection. 

An extension to the original chameleon model is to 
include an interaction term between the chameleon and 
electromagnetic (EM) field of the form jj-^F^F^ , first 
suggested by Brax et al. 2]. This leads to interconver- 
sion between chameleons and photons in the presence of 
a magnetic field, which alters the intensity and polariza- 



*Electronic address: a.c.davis@damtp.cam.ac.uk 
^Electronic address: C.A.0.Schelpe@damtp.cam.ac.uk 
t Electronic address: D.Shaw@qmul.ac.uk 



tion of any radiation entering the magnetic region. A 
coupling between the chameleon and EM fields gives rise 
to observational effects that could potentially be detected 
either astrophysically or in the laboratory @,i,l3. These 
predictions have lead to constraints on the mass parame- 
ter Mp describing the coupling strength to the EM field: 
M F > 2 x 10 6 GeV and Mp > 1.1 x 10 9 GeV Q. 

It was noted in Ref. Q that in low density back- 
grounds, such as the interstellar and intra-cluster medi- 
ums, the mixing between the chameleon field and pho- 
tons is indistinguishable from the mixing of any very 
light scalar field with a linear coupling to F^F^ U or 
a very light pseudo-scalar with a linear coupling to 
tjivpoF^ ' F pa . Such fields are collectively referred to as 
axion-like particles (ALPs). Standard ALPs have the 
same mass and photon-coupling everywhere. The best 
constraints on their parameters come from limits on their 
production in relatively high density regions such as stel- 
lar cores or in supernova 1987A [h\. Chameleon scalar 
fields represent one example of a subclass of ALPs that 
we term chameleonic ALPs. The defining feature of 
chameleonic ALPs is that the astrophysical constraints 
from high density regions are evaded or exponentially 
weakened. This may be because as the ambient density 
increases the field mass rises as in the chameleon model, 
or because the effective photon coupling weakens d, [f| . 
This latter possibility is realized in the Olive-Pospelov 
model for a density-dependent fine-structure constant 0] • 
The results we derive in this paper technically apply to 
all ALPs, however the parameter range for which they 
apply is strongly ruled out for all but chameleonic ALPs. 
The model for chameleon-like theories is described in de- 
tail in subsection II Al below. 
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In this paper we consider the propagation of the Cos- 
mic Microwave Background (CMB) radiation through 
the magnetic field of a galaxy cluster. We predict the 
modification to the CMB intensity and polarization in 
the direction of galaxy clusters arising from mixing with 
a chameleonic ALP. 

One of the significant foreground effects acting on the 
CMB as it passes through a galaxy cluster, which has 
been extensively verified by experiment, is the Sunyaev- 
Zel'dovich (SZ) effect. This arises from scattering of 
CMB photons off electrons in the galaxy cluster atmo- 
sphere, causing a change to the CMB intensity in the di- 
rection of the cluster. Measurements of the SZ effect can 
be compared to predictions of the CMB intensity decre- 
ment arising from the chameleon, and thus constrain the 
chameleon model parameters. 

This paper is organized as follows: in £ fl"fl the con- 
version between photons and a general scalar ALP in a 
magnetic field is analysed, and note how it depends on 
the structure of the magnetic field. In mill wc briefly 
review the observational evidence for magnetic fields in 
the intra-cluster medium (ICM) and detail their typical 
properties. In mVl we describe the standard SZ effect and 
show how this is modified when photons are converted 
to (chameleonic) ALPs in the intra-cluster medium. We 
note that the frequency and electron density dependence 
of the chameleonic SZ effect differs from that of the ther- 
mal and other standard contributions to the SZ effect. 
In fJV] we exploit this difference to constrain the prob- 
ability of photon-scalar conversion in the ICM. Under 
certain reasonable assumptions about the magnetic field 
structure, we convert these constraints to bounds on the 
effective photon to scalar coupling. In Wll we discuss the 
modification to the CMB polarization from chameleon- 
mixing. We summarize our results in Will 



although we do not require that B m = Bp. Varying this 
action with respect to tj> gives the scalar field equation: 



D6 = 



where 



cffl 



F z T 



dV cS . 



v{4>) 



■ F 2 T ) 



(2) 



B m (4>)T n 



Here T m = p m — 3P m , is the trace of the energy- 
momentum tensor and corresponds to the physical, mea- 
sured density or pressure. For this derivation we must 
remember that the particle masses are constant and in- 
dependent of <fi in the conformal Jordan frame with met- 
ric g^l = Bfg^y, while the length scales are measured in 
the Einstein conformal frame. Thus p nl = Bf p Jordan- 

Whether or not a given scalar-tensor theory is 
chameleon-like is determined by the form of its self- 
interaction potential, V(4>), and its coupling functions, 
B m (cf>/M) and Bp((f)/M). A simple and very general pre- 
scription is that the magnitude of small perturbations in 
In Bp(cj)/M) and In Bi(4>/M) about their values at the 
minimum of V c s{4>), and sourced by some small pertur- 
bation in T m or F 2 , decreases as T m becomes large and 
positive, i.e. as the ambient density increases. This im- 
plies that for any combination of i,j = F, m, 



Hij(F 2 ,T m ) 



(In Bi) t4 , (0 ml „/M) (In Bj ) i4> (r/w/M) 



V 



eff , 



(3) 



is a decreasing function of T m , where </Wn is given by 
V°6 (0mi n ; F 2 , T m ) = 0. We must also require that </> m i n 



exists and that VS5( 



;F 2 ,T m ) > for stability. 



A. Chameleon-like Theories 

Chameleon-like-theories, which include the standard 
chameleon model introduced by Khoury and Weltman 
[l| as well as the Olive-Pospelov, varying a cm model 0, 
are essentially generalized scalar-tensor theories and as 
such are described by the following action: 

r / m 2 i 

S = J d^x^—g (-^-K - -g^d^dA (1) 

-V{4>) - ^B F {<t>/M)F llv F^ 

where the ipi are the matter fields and <S ma tter is the 
matter action (excluding the kinetic term of electromag- 
netism). The Bi(<fi/M) determine the coupling of the 
scalar, (f>, to different matter species; Bp((j)/M) deter- 
mines the photon-scalar coupling. For simplicity we as- 
sume a universal matter coupling Bi(<f)/M) — B m (4>/M), 



1. The Chameleon Model 

In a standard chameleon theory, we assume that 
<p/M <C 1 and approximate Bp w 1 + c/)/M c g and 
B 2 ^ w 1 + 2(/>/M, where this defines M and M eS . Wc 
expect M er f ~ O(M) but this does not necessarily have 
to be the case. Assuming T m /M > |F 2 |/4A/ off , the con- 
dition on Hij is then equivalent to 

V(j>4>4> < 0. 

This implies that the (j> field equations will be non-linear 
in <f>. Additionally the existence of </>rni n and stability 
implies 

< 0, V M > 0. 

A typical choice of potential can always be described (in 
the limit of (f> > A) by 

. A"+ 4 
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where for reasons of naturalness A ~ O(Aq). For the 
chameleon field to be a suitable candidate for dark en- 
ergy, 

A = (2.4 ±0.3) x l(T 3 eV. 

A corollary of Eq. © in chameleon theories is that the 
chameleon field is heavier in dense environments than it 
is in sparse ones. From Eq. ([2]), we can see that the 
mass of small perturbations, m^, in the chameleon field 
is simply given by mi = V^. The chameleon field, <j>, 
rolls along the effective potential, V g, and the shape of 
V e f{ depends on the ambient density through T m . For 
baryonic matter under normal conditions, T m sw p m , and 
so when, as is often the case, p m ^> \F 2 \ and <f>/M -C 1 
we have, 



lei 



V{cj>) 



1 



M 



A plot of this effective potential, for when V{4>) is given 
by Eq. ([J) and n — 1, is shown in Fig. [TJ From this we 



v eff (o, P) 




^^p(1+<|>/M) 



FIG. 1: The chameleon effective potential, V e s, is the sum of 
the scalar potential, V(<f>), and a density-dependent term [lj. 

can see that the value of 4> at the minimum will decrease 
as the density of the surrounding matter, p m , increases, 
while the curvature of the effective potential at the min- 
imum will increase. It is clear from this that the mass 
of the chameleon is greater in relatively high density en- 
vironments such as laboratories on Earth, while in space 
the chameleon field is much lighter. It is this property 
that allows the chameleon field to have a strong coupling 
to matter, yet be light enough in space to produce in- 
teresting phenomena, such as non-negligible mixing with 
photons, whilst being heavy enough in the laboratory to 
avoid being ruled out by constraints on fifth-forces. A 
detailed discussion of the allowed chameleon parameters 
is given in Ref. [8|. 

The best direct upper-bounds on the matter cou- 
pling l/M come simply from the requirement that the 
chameleon field makes a negligible contribution to stan- 
dard quantum amplitudes, roughly l/M < 10~ 4 GeV -1 
[1, §. Stronger constraints have been derived on the 



photon coupling g e g — 1/M e g, Laser-based labora- 
tory searches for general ALPs, such as PVLAS [nj, 
and chameleons, such as GammeV [llj . give g e g < 
10 -6 GeV . Constraints on the production of starlight 
polarization from the mixing of photons with chameleon- 
like particles in the galactic magnetic field currently pro- 
vide the strongest upper-bound of g c g < 9 x 1CP 10 GeV -1 
[J] . With such couplings and for a standard choice of po- 
tential such as Eq. g| with n ~ 0(1), < KT 14 eV 
in the galactic and intra-cluster media. 



2. The Olive-Pospelov Model 

In any theory with the action of Eq. |T]), the ef- 
fective fine-structure constant, a cm , is proportional to 
1/ 1 Bf{4>/M). Observations of QSO spectra suggest that 
the a em in dust clouds roughly lOGyrs ago differs by a 
few parts in 10 6 from its value today as measured in the 
laborator y 11211 . This could be explained by either a time 
variation [l3|| in a cm , or a density dependence [?] since 
the clouds are much sparser than the laboratory environ- 
ment. In standard chameleon theories the variation in 
a em is constrained to be well below the level that can be 
detected, and so cannot explain this observation. Olive 
and Pospelov suggested a chameleon-like model in which 
a 0(1O -6 ) change in a cm (or any other constant of na- 
ture) from dense to very sparse environments is feasible. 
This model has, 



B F = 1 



B m = 1 



M 



2 V M 



where £p is an 0(1) constant. Whilst one may con- 
sider more complicated choices for the potential, a simple 
choice satisfying all the required conditions is V(<p) — 
const + m?(j) 2 j 2. In low density environments such as 
the interstellar or intra-cluster media m i n ~ 0, while 
in high density environments such as the laboratory 
4>ram ~ ^m- The total fractional change in a om from 
high to low density environments in this model is then 
8a/a ~ i F <fij2M 2 . 

The effective photon coupling in low density environ- 
ments is 



g cS = l/M eS ={hxB F ) A {(j> miu /M) 



10 



M 



10" 



1/2 



In high density backgrounds g c e is much smaller. In 
low density backgrounds, the mass of the Olive-Pospelov 
scalar, ~ m 2 and the constraints found in Ref. Q 
imply ?7i0 <C 10~ 14 eV. For the Olive-Pospelov model to 
satisfy all current experimental constraints, l/M must 
be neither too large nor too small. Olive and Pospelov 
gave the lower-bound l/M ^> 10 _9 GeV _1 , which for 
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8a/ a ~ 0(1(T 6 ) implies g eS > 10" 12 GeV -1 . Labora- 
tory tests such as GammeV and PVLAS are insensitive 
to OP scalar fields, because at their relatively high oper- 
ating densities the scalar-photon coupling is greatly sup- 
pressed. The best upper-bound on g c ff comes, as with 
standard chameleon fields, from limits on the produc- 
tion of starlight polarization in the galaxy. Thus if the 
OP model is to explain a 1CP 6 change in a om as sug- 
gested by observations of QSO spectra, we must have 
1(T 12 GeV- 1 « 9cff <9x 1CT 10 GeV. 



II. OPTICS WITH A SCALAR AXION-LIKE 
PARTICLE 

In the presence of an axion-like particle (ALP), 
the properties of electromagnetic radiation propagating 
through a magnetic field are altered as a result of mixing 
between the ALP and photons. Over the past 25 years or 
so, there has been agreat deal of work on this subject; 
for example Refs. [J [3, EE HI] or see Ref. [j| for a 
recent review. 

In this section, we present the equations which describe 
the mixing of a chameleon-like particle with the photon 
field in the presence of a magnetic field. We also present 
the solution of these equations, in the limit where the 
mixing is weak, for two different models of the magnetic 
field structure in the ICM. 

Varying the action in Eq. (Q} with respect to both <f> 
and Au, we find Eq. @ and 

v P [B F (d)/M)F^} = r, 

where is the background electromagnetic 4-current; 
V M J" - 0. 

We are concerned with the propagation of light in 
sparse astrophysical backgrounds containing a magnetic 
field B(a; M ). The background values of the chameleon 
and photon fields are assumed to be slowly varying over 
length and time scales of 0(1/uj) where ui is the proper 
frequency of the electromagnetic radiation being consid- 
ered. We define the perturbations in the photon and 
chameleon fields about their background values to be, 
fl/j = — A^ and ip = (ft — (ft, respectively. When elec- 
tromagnetic radiation propagates through a plasma with 
electron number density n Cl an effective photon mass- 
squared of = 4Tra om n e /m e is induced; w p i is the 
plasma frequency. Following Ref. [H in ignoring terms 
that are second order in the perturbations ip and a M , and 
assuming uj 3> H where H is the Hubble parameter, we 
find 

_ 2 Vip x B 2 
a+V 2 a ~ — b h^,a, 



-<p + v> 



M eff 

B • (V x a) 



cir 



where 1/M e g = B F ^((ft) and m| = V^s^^cf), T m , A^). 



We consider a photon field propagating in the z direc- 
tion, so that in an orthonormal Cartesian basis (x, y, z) 
we have a = (7^, j y , 0) T . The above equations of mo- 
tion for the chameleon and photon polarization states 
can then be written in matrix form: 



-pi 




d z 





M.rt Uz 



M ca °z 1 






( lx 










ml J 









= 0. 



Following a similar procedure to that in Ref. [l5| . 
we assume the fields vary slowly over time and that 
the refractive index is close to unity which requires 
m 2 ,/2w 2 , u> 2 1 /2uj 2 and \B\/2wM eS all < 1. Defining 

7, = ^ z ) e Mz~t)+i(i(z) and ^ = p e iu(z-t)+i0(z) where 

/3 Z = — cjp[(z)/2w, we approximate — d 2 « cu 2 and 
uj 2 + d 2 « 2u>(uj + id z ). In addition to this, we define 
the state vector 



u 



e 2tA ^p(z) 



where 



A(z) 



4lo 



-dx. 



and m 2 ff = m 2 ^{z) — cl> 2 [(z), which simplifies the above 
mixing field equations for a and (p to 



B(z) 



2M, 



where 



B(z) 






B y e 2lA 



fir 






-B r e 2iA 



(5) 



-B y e- 2tA 
B x e~ 2lA 




To solve this system of equations we expand, u = Uo 
Ui + u 2 + . . ., such that 



B[z) 



i+i 



2M, 



off 



and neglect the higher order terms from mixing. We 
define 



Mi = 



B{x) 



2M, 



dx. 



eff 



We say that mixing is weak when 

tr (m\(z)Mi(z)) < 1. 

Thus Eq. ((5|) can be solved in the limit of weak-mixing: 
u(z) ~ [I + Mi(z) + M 2 (z)} u(0), 
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where M2 = Jq M' 1 (x)Mi(x)dx, which can be written In what follows we assume that the scalar field is very 



u(z) 



l + M 1 (z) + -(Mc(z)+Ml(z)) 



u(0), 



where explicitly 



and 



Mi 



Mr 



Ai = 



— A* 
A* 




Ay Aj 

— C _r<* 
^yy ^'xy 

Cxy C X x 





Cy 



dx 



Bi(x)e 



2iA(x) 



2M 



cir 



dx (A*'(x)Aj(x) - A*(x)A' j (x)) 



We note that M\ = —M n - 

The polarization state of radiation is described by 
its Stokes parameters: intensity I- 1 {z) — \^ix{z)\ 2 + 
\j y (z)\ 2 , linear polarization Q(z) = I'Jxiz)] 2 — \ly( z )\ 2 
and U(z) = 2He (j*(z)"f y (z)), and circular polarization 
V(z ) = 21m (7* (z)jy(z)). Assuming there is no initial 
chameleon flux, Ia = \lp\ 2 = 0, then to leading order the 
final photon intensity is given by 

L ( (z) = J 7 (0)(1-7V*(*)) 

+Q(0)Q q (z) + U{0)Q u {z) + V{0)Q y (z), 

where we have defined, 



Qu(z) 
Qy{z) 



(\A x {z)\ 2 + \A y {zt 
(\A x (z)\ 2 -\A v (zf 



MA* x A y ), 
Im{A*A y ). 



light so that m 2 s 



Spatial variations in 



m 2 s are then entirely due to spatial variations in n c . 

The magnetic fields of galaxies and galaxy clus- 
ters have been measured to have a regular component 
B reg (z), with a typical length scale of variation simi- 
lar to the size of the object, and a turbulent component 
(5B(z) which undergoes 0(1) variations and reversals on 
much smaller scales. It is common, both in the literature 
concerning photon- ALP conversion in galaxy and cluster 
magnetic fields and in that concerning the measurement 
of those fields, to use a simple cell model to describe <5B. 
In this cell model one defines L C oh to be the length scale 
over which the random magnetic field is coherent, i.e. 
the typical scale over which field reversals of <5B occur. 
A given path of length L is then divided into N magnetic 
domains, N = L/L co h. The magnetic field strength and 
electron density are assumed to be constant along the 
path, but the random magnetic field component, ^B x , is 
taken to have a different random orientation in each of 
the iV magnetic domains. Although the cell model for <5B 
is common in the literature and reveals most of the key 
features associated with astrophysical photon- ALP mix- 
ing, it is not accurate especially when |A(z = L)\ 3> 1. 
Additionally, as was first noted in Ref. [17| , in this model 
we do not have V • #B = as we should. A more real- 
istic model for B is to assume a spectrum of fluctua- 
tions running from some very large scale (e.g. the scale 
of the galaxy or cluster) down to some very small scale 
(<C i C oh)- We describe these fluctuations by a correlation 
function i?s(x) = (<5B(y)5B(x + y)) and the associated 
power spectrum Pe(fc); here the angled brackets indicate 
the expectation of the quantity inside them. We can, and 
should, also allow for a similar spectrum of fluctuations 
in the electron number density with some power spec- 
trum P^{k). We now evaluate V 1 ^ ( j > and the Qi in both 
this power spectrum model and the cell model. 



In this article we consider mixing with the CMB, 
the intrinsic polarization of which is small, i.e. 
v/Q 2 (0) + U 2 (0) + V 2 (0)/I 7 {0) < 1. It follows that to 
leading order the modification to the Stokes parameters 
of the CMB radiation, from the effect of photon-ALP 



mixing, is 



AS(z) = A 




-Vj^iz) 

Qq(*) 
Qu{z) 

QAz) 



J 7 (0). 



To evaluate ? 7 „^ and the Qi we need a model for the 
spatial variation of the magnetic field and effective mass. 
Typically n c ~ 10 -3 -10~ 2 cm~ 3 in galaxies and galaxy 



clusters and so uj p \ ~ 0(10 ) eV. 
like-theories detailed in £jl A[ it follows that cu 2 , ^S> m 2 



For the chameleon- 
pi 



The Cell Model 



In the cell model, we take m 2 s = const and B = B reg + 
<5B where B rcg ps const over the path length L. For <5B 
we take (<5B)a; = -Brand cos 9 n and (SJi) y = i? ran d sin # n 
for (n— l)L C oh < z < rtL con where the n are independent 
random variables with identical distributions £7(0, 27r]. 
We then have, over a path length L, 



A,. = 



A — 



(Areg), + (5A)i, 

2B rc ^e lA 



■ sin A, 



M cS m 2 



eff 



N-l 

E 

r=0 



n l >e N 
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where A = m 2 s L/4:Uj and = (cos 9 n , sin#„, 0) T . It 
follows that the expected values of P 7< _,0 and the Qi are, 



p. 



(p 7 



1 / 2P rcg tj 

2 ^M off m2 ff 



sin 2 A (6) 



JV / 2B rand w 



o , sin . . . 
2 / \ N 



2 \Meffm* ff 



2 \ M eS m 2 s 



2 V M cfrm 2 ff 



Qv = 0, 



sity fluctuations, P B (k) and P^(k), are defined: 
I 



47T 



/ d 3 fce 2 ™ kx 



k 2 dkP B (k) 



PB(k) 
sin(27rfca;) 



2irkx 

R N ( X ) = J- f d 3 fce 2 ™ kx P N (fc) 

47T J 

sin(2irkx) 



k z dkP N (k)- 



2nkx 



Following Ref. [18| . we define correlation length scales, 
Lb and Ln, for the magnetic and electron density fluc- 
tuations respectively: 



where we have used (<5B) = 0, and defined (B rcg ) a; = 
.Breg cos # reg and (B reg )j/ = B reg sin# reg . We also note 
that for A/N > 1 and TV > I, the variance of the Qi 
at a given frequency are 0(P 2 ^). We present the de- 
tailed evaluation of P 7< _>0 and the Qi in the cell model in 
Appendix IA II These results were first presented in Ref. 



Apt 



We note that if, as is generally the case, Brand ~ 
0(P rog ), the effect of multiple box-like magnetic regions 
with coherence length L corl <C L is to enhance the ex- 
pected photon- ALP mixing probability by a factor of N. 



B. The Power Spectrum Model 

A more generalised prescription for describing the mag- 
netic field and electron density fluctuations is with a 
power spectrum model. As before, m 2 s ~ ~ UJ p\ x n <=- 
We assume that, for each i and for fixed x, the turbulent 
component of the magnetic field #B j (x) is approximately 
a Gaussian random variable. We assume that the fluctu- 
ations are isotropic and so 



i? By -(x;y) = (SB^y^B^x + y)) = ^R B (ic;y)6 ij . 



Finally we require (at least approximate) position inde- 
pendence for the fluctuations so that B B (x; y) ~ R B (x). 
Rb(x) is the auto-correlation function for the magnetic 
field fluctuations. The electron number density is divided 
into a constant part and a fluctuating part, n c = h c +Sn e , 
and we assume that 1 + 5n c /h c is a log-normally dis- 
tributed random variable with mean 1 and variance (<5 2 ), 
where S n = 8n c /n c . We define the electron density auto- 
correlation function by, 

Rn{x) = (Sn c (y)Sn c (x + y)) . 
The power spectra for the magnetic and electron den- 



Lb/n — 



/ °°fcdfcP B/ N(fc) 

2/ ° o FdfcP B/N (fc)- 



(7) 



We define A = m 2 s L/Auj ~ —ATra em n e L/4m e u>, and 



|A/ttL|. In Appendix lA~2l we evaluate P. 



at 



low frequencies, i.e. large |A|, in the power spectrum 



model. We assume there is some k 



(B/N) 



<C fecrit Such 



that Pb/n(^) decreases faster than k as k — ► oo for 
k > k^-l^ ■ This is equivalent to assuming that the dom- 
inant contribution to (<5B 2 ) and (n 2 ) comes from spatial 
scales than are much larger than fc~ it . We expect that 
this dominant contribution will come from scales of the 
order of the coherence lengths of the magnetic field and 
electron density fluctuations and so are assuming that 
k~ lit -C L B , Pnj. This and the other assumptions given 
above lead to, 



(P 7 



b) ~ 2 

2P reg w 
M cS m 2 cS 



2B c ffUi 



M cS m* s 



r 5 



;(2A) 



D 



2 T 



L 



2AM 2 

err 



(8) 



4w B (k CIit ). 



where 



J N = l + (S 2 n ), 



W B/N (k crit ) 



^((zxB) 2 ) 
P 2 cg /2 + (SB 2 ) /3, 
kdk P B / N (fc). 

fccrit 



It is straightforward to check that Q v = 0, and 

Q q ~ Qocos20 reg , 
Qu ~ Qosin20 reg , 
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where 



Qo 



1 / 2B rcE uj 
4 \M eS m 2 cS 

_1 / 2P reg o; 
4 ^M eff m 2 ff 



r 3 



;(2A) 



B icg L 

16M 2 , 



efif"e 

We note that, in contrast to the cell model, here we have 
an extra term in the expressions for the mean induced 
Stokes parameters which scales as L. This term is asso- 
ciated with electron density fluctuations and is usually 
dominant when A ^> 1 , as is the case for CMB photons 
in galaxy clusters. 

We have assumed that near k = fc cr it and for all large 
fc, both Pe(fc) and Pn(&) are decreasing faster than fc~ 3 . 
We further assume that nCeir k — ^crit i Pe(fc) oc P^(k) oc 
f or some a < _i define normalization con- 

stants Ck and Cn: 



k 2 P B (k) = 2C K 



k 



(9) 



fc 2 P N (fc) = 2(27r) 1 / 3 C*Sfc Q , (10) 

where fco = lkpc -1 . These forms hold for k ~ 0(fc cr it)- 
We have chosen the normalization for Pb and Pn to 
be consistent with the conventions of Refs. and 
[20l ] respectively, where a = —5/3 corresponds to three- 
dimensional Kolmogorov turbulence. It follows that, 



W B (fc cr it) 
W N (fccrit) 

We have 

A£i«2.4x 10- 2 P c 



2Ck / fceri 



|a| \ fc 
2(2^)!/3 



100 GHz 



°N K crit- 



10 3 cm 



> (11) 



where v = uj/2-k is the frequency of the electromagnetic 
radiation. For galaxy clusters n c ~ 10 _3 -10 -2 cm~ 3 
and for CMB photons v ~ 30 - 300 GHz. The typical 
critical length scale is then much smaller than a parsec, 
k~ r l t w 10 _3 -0.1pc. The form and magnitude of Pe(fc) 
and P]yr(fc) in galaxy clusters for such small scales are not 
known empirically. 

We note that the cell model is in fact a sub-class of 
the power spectrum model with an appropriate choice 
for Wb and W^. When referring to the power spectrum 
model in what follows we will be referring to the specific 
case of Pe(fc) oc Psr(fc) oc k a ~ 2 at large fc. 



III. MAGNETIC FIELDS IN GALAXY 
CLUSTERS 



of the magnetic field along the line of sight. As the CMB 
photons propagate towards Earth they pass through clus- 
ters of galaxies. The space between individual galaxies 
in the cluster is filled with a hot plasma which supports a 
magnetic field. In this article, we are focused on the mix- 
ing of CMB photons and light (chameleonic) scalars in 
the intra-cluster (IC) medium. In this section we briefly 
review the observed properties of IC magnetic fields. 

Intra-cluster magnetic fields have been observed to 
have field strengths as high as 30 fiG [2l[ . The galaxy 
clusters can be divided into two types: non-cooling core 
clusters for which the magnetic field in the ICM is gener- 
ally found to be a few to 10 /xG and ordered on scales of 
10 — 20 kpc; while in cooling core systems the magnetic 
field strength is generally higher (around 10 — 40 /iG) and 
ordered on smaller scales of a few to 10 kpc [22|. In Ta- 
ble U we list the typical ranges for some of the galaxy 
cluster properties. 

One of the common methods for determining cluster 
magnetic fields is the observation of Faraday rotation 
measures (RMs). The RM corresponds to a rotation of 
the polarization angle of an intrinsically polarized radio 
source located behind or embedded within the galaxy 
cluster. The RM along a line of sight in the z direction 
of a source at z = z s is given by: 

f Zs 

RM(z s z) = ao / n 6 (xz)B\\ [x%)Ax, 
Jo 

where ao = a^/ir^-^m 2 , B\\ is the magnetic field along 
the line of sight, and the observer is at x = 0. For a 
uniform magnetic field over a length L along the line 
of sight and in a region of constant density, we expect 
(RM) = cioB\\n c L. However, the dispersion of Faraday 
RMs from extended radio sources suggests that this is 
not a realistic model for IC magnetic fields. Both the 
dispersion of the measured RMs and simulations suggest 
that IC magnetic fields typically have a sizeable compo- 
nent that is tangled on scales smaller than the cluster. 
The simplest model for a tangled magnetic field is, as we 
mentioned in the previous section, the cell model where 
the cluster is divided up into uniform cells with sides of 
length i C oh- The total magnetic field strength is assumed 
to be the same in every cell, but the orientation of the 
magnetic field in each cell is random. In the cell model 
the RM is given by a Gaussian distribution with zero 
mean, and variance 



(RM 2 ) 



'Z T 

u c °h 



(B 2 (z))n 2 e (z)dz, 



where (B 2 ) is the average value of B 2 along the line 
of sight. We noted in the previous section that a more 
realistic model for the magnetic field is to assume a power 
spectrum of fluctuations. In this model one finds |18| . 



We saw in the previous section that the mixing of pho- 
tons with light scalars depends crucially on the structure 



(RM 2 ) 



(B 2 (z))n 2 (^)dz. 
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Features of a typical galaxy cluster : 




magnetic field strength, (B 2 ) 1 ^ 2 
coherence length of magnetic field, L co h 
number of magnetic domains, N = L/L co h 
electron number density, n c 


0.1 - 30/^0, 
1 - 100 kpc, 
100 - 10 3 

10" 4 - 10" 2 cm -3 , 



TABLE I: List of the parameters for typical galaxy clusters [23l. |24||. 



To match the same observations with the same total held 
strength, Lb in the power spectrum model should be 
equal to L co h in the cell model. Whichever model is used, 
determination of the magnetic field from Faraday rota- 
tion measures is heavily reliant on the estimated coher- 
ence length of the magnetic field. A study by Clarke et 
al. (2f| examining the Faraday RMs of 16 "normal" clus- 
ters, i.e. free of strong radio halos and widespread cooling 
flows, found the following relationship between the aver- 
age field strength and the inferred coherence length for 
the IC magnetic field: 

(l B D^ = fi - 10 (^^) V2 ^VG (12) 

where /175 is defined by Hq — 75 kms _1 Mpc _1 . This 
study excludes the cooling core systems which generally 
possess slightly higher magnetic field strengths. 

The ambient electron number density also plays a cru- 
cial role in determining the strength of photon-scalar 
mixing. The electron density of cluster atmospheres can 
be determined from X-ray surface brightness observa- 
tions, for example using ROSAT, which are fitted to a 
/3 profile (see Ref. (H): 

/ r 2\ -3/3/2 

n e (r) = n ^1 + -^J , 

where /3 ~ O(l) and positive. The magnetic field is ex- 
pected to decline with radius in a similar fashion to the 
electron density. Magneto-hydrodynamic (MHD) simula- 

1/2 

tions of galaxy cluster formation suggest that (B 2 ) cx 
(n e } v , 77 = 1 while observations of Abell 119 suggest 
77 = 0.9 [l?) • Two simple theoretical arguments for the 
value of 77, are (1) that the magnetic field is frozen-in 
and hence 77 = 2/3, or (2) the total magnetic energy, 
cx (B 2 ), scales as the thermal energy, cx n e , and hence 
77 = 1/2 (26j. Faraday RM observations of clusters (25| 
have found that the magnetic field extends as far as the 
ROSAT-detectable X-ray emission (coming purely from 
the electron density distribution). Typically this corre- 
sponds to a scale of ~ 500 kpc. 

A. Power Spectra of Cluster Magnetic Fields 

As discussed in miBl fluctuations across a range of 
scales in the magnetic field and electron density in clus- 
ters can be described by their respective power spectra, 



PB(k) and P^(k). On small scales the power spectra 
are parameterised by a power law form as given by Eqs. 
© and (fTU]) . In these equations the index a = —5/3 
corresponds to the special case of three dimensional Kol- 
mogorov turbulence. According to Kolmogorov's 1941 
theory of turbulence, the power spectrum of three di- 
mensional turbulence on the smallest scales is universal 
and has P(k) cx fc -11 / 3 . 

The values of the normalization constants, Cj^ and 
Ck, in our own galaxy have been determined by Arm- 
strong et al. and Minter & Spangler [28[ respec- 
tively. Armstrong et al. found a ~ —5/3 for spatial 
scales 10 6 m < k~ x < 10 13 m, and C§ w lO^m" 20 / 3 . 
Minter & Spangler also found a w —5/3 for spatial 
scales 0.01 pc < fc" 1 < 3.6 pc, and C K w 0.95 x 
10~ 12 erg cm~ 3 kpc. On larger scales Han et al. [l9[ found 
that the magnetic fluctuation power spectrum flattens 
with a « -0.37 ± 0.1 for 0.5 kpc < fc" 1 < 15 kpc, and 
Ck ~ 6.8 x 10 -13 ergcm~ 3 kpc on these scales. 

Enfilin & Vogt [17| and Murgia et al. [18| showed that 
information about the power spectrum of galaxy cluster 
magnetic fluctuations could also be gained from detailed 
RM images. The inferred galaxy cluster magnetic power 
spectra are consistent with a power law form and — 2 < 
a < in different regions of the clusters. For instance, 
in Ref. [29( Vogt & Enfilin analysed Faraday rotation 
maps of three galaxy clusters: Abell 400, Abell 2634 and 
Hydra A. They found that for k > lkpc -1 , the power 
spectra of all three clusters were consistent with a power 
law form and —2 < a < —1.6, and hence consistent with 
a Kolomogorov power spectrum (a — —5/3) on small 
scales. For k < lkpc -1 , fc 2 Pe was found to be roughly 
flat or slightly increasing with fc. 

Such measurements only probe the power spectrum at 
spatial scales larger than a few kiloparsecs. In this work, 
we find that the photon to scalar field conversion rate is 
sensitive to Pb(^) and Pu(k) on scales ~ fc"^ (Eq. (fTTj) h 
Typically k~] t « 10 _3 -0.1 pc for the CMB in galaxy clus- 
ters. There are no observations of the form or magnitude 
of -Pb(^) and ffo(fc) on such small scales. We therefore es- 
timate a by assuming that on scales < fc"^ , we have three 
dimensional Kolmogorov turbulence. The dominant con- 
tribution to (<5B 2 ) comes from scales where Pb switches 
from dropping off more slowly than fc -3 to dropping off 
faster than fc~ 3 as k increases. We define this scale to be 
fc* . To estimate Ck and Cn we assume a simple model for 
P b /n(^) whereby for fc > fc* up to and including fc cr it, we 
have Kolmogorov turbulence, and for fc < fc*, the power 
spectrum decreases more slowly than fc~ 3 with increasing 
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k. We also assume that A; 2 Pb/n(^) — ¥ const as k — > so 
that the coherence lengths L B and Ln (Eq. 0) are well 
defined. We therefore assume, 

-5/3 



fc 2 P B (fc) « 2C K ( — ) , k>K 



in line with observations. Following a similar procedure 



for we estimate, 



2{2ir) 1/3 C^n- 2 w (0.27-0.45) 



<K 2 > 



-2/3 



where (fc 2 ) /n 2 = Jn — 1, and we approximate Ln by 
Lb- The above expressions provide order of magnitude 
estimates for Ck and 



where /3 > — 1, and we can ignore the contribution of 
k 2 PB(k) for k < fc m i n . It follows that, 



<<5B 2 ) 



1 



d 3 kP B (k) 



47T 

[2/i-hs(aO+3]C K fc, 



kn 



-5/3 



i B (<5B 2 ) = i J kdkP B (k) 



' 1 



-5/3 



where x = fc m i n /fc* and we have defined 



1 



From which we find, 



and 



where 



k _ L -i i ftfcHj 



kTc K = ^ m h(x){SB 2 )L B 2/3 , 



h{x) 



tf'%fp{x) + |) 2 / 3 
4(2/ 1+(3 ( a; )+ 3)5/3- 



On large scales, 100 k pc" 1 < k < lOkpc -1 , observations 
such as those in Ref. [29| are consistent with /3 ss 0. We 
approximate /c m i n ~ ^(-^a us t) where L c i ust is the length 
scale of the magnetic region and we have fc* ~ 0{L^ 1 ). 
Thus k*/k m i n ~ 10 — 200 which gives, taking /3 ~ 0, 



Hence, 



■ hl^/fcmin) + 



2/3 



1.5-2.6. 



-2/3 



2fcg /J C K S3 (0.27-0.45) (SB 2 ) L 

For the galactic magnetic field where (<5B 2 ) 
and Lb ~ 20-50 pc, this gives 



1/2 



3 /LtG, 



C K ~ (0.7-2.2) x 10~ 12 ergcuT 3 kpc, 



IV. CHAMELEONIC CORRECTIONS TO THE 
SUNYAEV-ZEL'DOVICH EFFECT 

One of the well-documented foreground effects acting 
on the CMB as it passes through galaxy clusters is the 
Sunyaev-Zel'dovich (SZ) effect. It arises from inverse- 
Compton scattering of CMB photons with free electrons 
in a hot plasma. The scattering process redistributes 
the energy of the photons causing a distortion in the fre- 
quency spectrum of the intensity. This is observed as 
a change in the apparent brightness of the CMB radia- 
tion in the direction of galaxy clusters. If a chameleon 
field or other light scalar couples to photons, the flux 
of CMB photons passing through galaxy clusters would 
additionally be altered by the mixing of this scalar field 
with CMB photons. As we shall see, for this new effect 
to be detectable the required effective scalar-photon cou- 
pling, g c ff, and effective scalar mass, m^, are ruled out 
for standard axion-like-particlcs. A standard ALP with 
the required properties would violate constraints on ALP 
production in the Sun, He burning stars and SN1987A. 
However, all of these constraints feature ALP produc- 
tion in high density regions. If the ALP is a chameleon- 
like scalar field, then constraints from high density re- 
gions only very weakly constrain the properties of the 
scalar field in low density regions such as the ICM. A 
chameleon-like scalar field includes the chameleon and 
OP models discussed earlier or any other non-standard 
ALP for which one or both of g e g and l/m^ decrease 
strongly enough as the ambient density increases. We 
refer to the correction to the Sunyaev-Zel'dovich effect 
due to mixing with some (chameleon-like) ALP as the 
chameleonic SZ (CSZ) effect. 

In this section, we briefly review the properties of the 
CMB and the standard thermal Sunyaev-Zel'dovich ef- 
fect before examining the chameleonic SZ effect. 



A. The CMB 

The Cosmic Microwave Background (CMB) was emit- 
ted in the early Universe from the surface of last scat- 
tering at the time of recombination. This radiation is 
approximately uniform across the sky and has the inten- 
sity spectrum of a black-body radiating at a temperature 
of T = 2.725A:: 

/cbTq uj 2 x 



In = 



2vr 2 



1 
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where x = uj/k^To, and fee is the Boltzmann constant. In 
addition to this monopolc there are small fluctuations in 
the primary CMB arising from primordial density fluc- 
tuations at the time of emission. Small fluctuations in 
intensity are related to the changes in the characteristic 
black-body temperature by 



6I_ 



5T 



T 



The root mean squared average of these primary fluctu- 
ations is 0(lO 2 )/iK. As the CMB radiation propagates 
towards Earth there are foreground effects modifying the 
intensity which give rise to secondary fluctuations in the 
CMB. We expect the chameleonic SZ effect in galaxy 
clusters to be one such foreground. When considering 
these foreground effects, one generally treats the primary 
CMB fluctuations as Gaussian noise on top of the CMB 
monopole signal. 

The primordial CMB radiation is also weakly polar- 
ized. The intensity and polarization of radiation in gen- 
eral is described by its Stokes vector (I, Q, U, V). The 
fractional linear polarization of the CMB is of the order 
(Q 2 ) 1 / 2 = (U 2 ) 1/2 ~ 10~ 6 , while the circular polariza- 
tion component, V , is predicted to be zero from parity 
considerations in the early Universe. 



B. The Thermal SZ Effect 

In the standard SZ effect the predicted fractional 
change to the temperature of CMB photons passing 
through a galaxy cluster is 30] 



ATI 



sz 



knT„ 



Tn 



-T() 



(xcoth(|)-4) , (13) 



where x — uj/k-gTo, T c is the temperature of the elec- 
trons in the plasma, and To = J ain e (l)dl is the optical 
depth given by the integral along the line of sight of the 
number density of electrons multiplied by the Thompson 
cross-section, or = 6.65 x 10~ 29 m 2 . At low frequencies 
the SZ effect is seen as a decrement in the expected tem- 
perature from the CMB, while at high frequencies it is 
seen as a boost in the temperature. There are relativistic 
corrections to the above formula, but for the specific case 
of the Coma cluster that we consider in this paper the 
peculiar velocity is small and higher order corrections to 
the thermal SZ effect are negligible [3lj . 



C. The Chameleonic SZ Effect 

We found in Sjn] that on passing through a magnetic 
field of length L, the photon-scalar mixing alters the in- 
tensity of light from Iq to If , where If is given by 



If = J (1-7V. 
+QoQ q (L) 



\-UoQ u (L) + VoQ v (0). 



The Qo, Uq and Vq are the initial values of the Stokes pa- 
rameters. The CMB radiation is only very weakly polar- 
ized and observations confirm (Qq) 1 ^ 2 /la, (Uq) 1 ^ 2 /Iq ~ 

O(10- 6 ) and (V 2 ) 1/2 /I a < 10~ 6 . The root mean 
squared values of P 7 „,(, and the Qi are, however, all of 
similar magnitude, 0{P 1 ^,j,). In addition to this the 
measurement issues arising from a finite frequency bin, 
as discussed in Appendix IA 3[ mean that the observed 
variance will be significantly reduced. Hence, for CMB 
photons, we have that 

If w/ (l-'P^(i)), 

and so the change in intensity due to the chameleonic SZ 
effect, A/csz, is given by 



A/csz 



Transforming this to a change in temperature we have, 



AT, 



csz 



1 



AI, 



csz 



1 



T 



Iq 



(L). 



where as before x = uj/kT . We evaluated the expected 
value of P 7< _,^(L) for different models of the IC magnetic 
field in fjll] above, and we note that it depends both on 
n Cl B and frequency u). We assume that inside the cluster 

1 /2 

(B 2 ) oc {n e ) n . This implies for the cell model we have, 



AI ( 



CSZ -2(l-7j) 2 



(14) 



whereas in the power spectrum model, where Pb(&) oc 
P N (k) cx k°- 2 for k > k crit (Eq. (JTTJ) ) and a < -1, we 
have to leading order, 



AI, 



(15) 



for the case a > —2. For a < —2, we find the above 
behaviour for higher frequencies, whereas at lower fre- 
quencies we find the same leading order scaling as we 
had in the cell model. A Kolmogorov spectrum for the 
magnetic turbulence has a — —5/3, if we take rj = 0.9 as 
suggested by observations, then the CSZ effect scales as 

2/15 5/3 

For comparison, the thermal SZ effect is proportional 
to n c and has a non-power law frequency dependence. 
We note that for a < —1 and f) < 1, A/csz/Io al- 
ways has a shallower n e dependence than A/sz/Io- in 
principle, the chameleonic SZ effect can be distinguished 
from the standard contributions to the SZ effect using 
both its electron density and frequency dependences. It 
should be noted, though, that in both Eqs. (fl"4"l) and (fl"5|) 
we have assumed that the magnetic and electron density 
correlation lengths do not vary as n c varies. Simulations 
suggest that the correlation length increases with radius, 
although the precise form of this variation is uncertain. If 
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we assume that L co h = Lb oc Ln, and estimate the high- 
k power spectrum according to the prescription given in 
gill Al we see that 



A7, 



csz 



c w ^coh > 



where ao = a if —2 < a < — 1, and «o = —2 in the cell 
model. The I/ co h scaling for a < — 2 is more complicated, 
at higher frequencies it is given by the above formula 
with ao = a, whereas at lower frequencies, Aicsz is 
independent of L C oh at leading order. In what follows 
we assume that — 2 < a < — 1 for the power spectrum 
model. We discuss how the spatial variation of A/csz 
can be used to constrain photon-scalar mixing further in 
WU below. 



V. APPLICATION 

In this section, we apply the results derived above 
to constrain the probability of photon-scalar mixing at 
CMB frequencies in the intra-cluster medium. We found 
that in the power spectrum model, 'P 7< -,^ oc tu~ a , for 
— 2 < a < — 1 where a is the slope of the power spec- 
tra for the magnetic and electron density fluctuations in 
the region k > k CY it (Eq. (|11|). The cell model for the 
magnetic field effectively has a — — 2. We write 



•,{L,u ) 



for some fixed frequency ojq of the order of the CMB 
frequency. For given objects we are able to constrain 
P 7 ^,p(L, u>q), the average probability of photon to scalar 
conversion at ujq. We pick uiq so that for the range, — 2 < 
a < — 1, the resulting constraint on V^^,^(L,iS) is almost 
independent of a.. 



A. Magnitude of Chameleonic SZ Effect 

Using our evaluation of "P~++$(L,iS) in Eqs. ([6]) and 
©, a constraint on P 7 «_>^(L, Wq) can be be transformed 
into a constraint on the effective photon-scalar coupling 
9cB = I/Meg. V-y^^iL^LOo) depends on quantities which 
are not known a priori, such as (B 2 ), (n c ), (u^), Lb, £n, 
Cb and Cn- The accuracy to which these parameters 
can be measured or estimated will therefore limit the 
constraints on g c g. The derived constraint on g e g will 
also depend strongly on a. We explicitly consider the 
relation between g g and V^^tf, for two models: the cell 
model where effectively a — —2 and the Kolmogorov 
turbulence model where a — —5/3. The constraint on 
V-y<^$(L, u>o) is, however, reasonably model independent, 
with only a slight dependence on a. 

In the cell model, Eq. © gives, 



V (ccll \L,is) 



(3.2 x 10 



10\ 



10/iG r 



0.01 



(-klkpcJ 



-1,,2 



1 /2 

where B W ^ G = (SB 2 ) /10/iG, ra .oi = n c /10" 2 cm" 3 , 
1/214 = 1^/214 GHz, L 200 - L/200kpc, £™kpc = 
L co h/lkpc and 510 = g c ff/(10~ 10 GeV -1 ). We have as- 
sumed A/TV ^ 1 and that measurements are over a finite 
frequency bin so that sin 2 (A//V) 1/2. 

A potentially more realistic model for the magnetic 
field structure is to assume that the magnetic and elec- 
tron density power spectra follow a Kolmogorov power 
law for k ranging from L~^ h to at least fc cr it, where fc cr it 
is given by Eq. (TTTj) . For k < L~ oh we assume that 
k 2 P{k) is cither approximately flat or increasing. The 
dominant contribution to the magnetic field fluctuations 
then comes from spatial scales of about L co h- In i jlll Al 
we considered a parameterisation of Pb(^) with these 
properties, and found 

2C K fco /3 = (0.27-0.45) (SB 2 ) L B 2/3 . 

We also assume that the Sn e /n e power spectrum is pro- 
portional to that of SB. Under these assumptions, Eq. 
© gives, 



7>( Kolm °)(L,z/) = (2.4-3.8 x 10" 7 ) g 2 w B 2 



10^G n 0.01 



-5/3 



5/3 
y 214 



/ \-2/3 ,j 

L200 [ L Tkpc) {- 



(2/a 



A reasonable value for Jn based on electron density 
fluctuations in our galaxy is In ~ 1-2, and hence 
I 2 (27 N - 1) ~ 1-12. 

The relationship between magnetic field strength and 
coherence length found from a study of 16 non-cooling 
core clusters (Eq. (TT2"])) allows us to simplify the above 
prediction for the case of non-cooling core clusters. For 
the range 1 kpc < L co h < 100 kpc, we find that in the cell 
model, 



.0 x 10" 14 < 



•p(cell) 



5l0^200« . 01^214 



< 3.2 x 10" 



Similarly in the Kolmogorov power spectrum model, in- 
cluding the range for In ~ 1-2, we find 



2.2 x 10" 



< 



■p(Kolmo) 



5/3 5/3 



< 3.8 x 10" 



Sio-^oo^o.oi v. 



Cooling core clusters generally have slightly higher 
magnetic fields and so we might expect them to exhibit a 
greater chameleon effect. Hydra A is a good example of 
a galaxy cluster with strong cooling flows [2 If] . Faraday 
RMs from an embedded radio source in the cluster have 
been analysed to suggest a smooth component to the 
magnetic field of strength 6 [iQ ordered on 100 kpc scales 
and a tangled component of strength 30 \xG ordered on 
scales of 4 kpc. It is this latter component which makes 
the dominant contribution to P 7 ^0. X-ray data of the 
cluster indicate a core radius of R corc = 130 kpc and aver- 
age electron density of 10~ 2 cm -3 . We approximate the 
path length, L, as being equal to twice the core radius 
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Experiment 


^/GHz 


Sp/GRz 


ATsz/fiK 


OVRO 


32.0 


6.5 


-520 ± 83 


WMAP 


60.8 


13.0 


-240 ± 180 


WMAP 


93.5 


19 


-340 ± 180 


MITO 


143 


30 


-184 ±39 


MITO 


214 


30 


-32 ± 79 


MITO 


272 


32 


172 ± 36 



SZ Effect for the Coma Cluster 



TABLE II: SZ measurements of the Coma cluster. 



Lrore = 2R cmc = 260 kpc. With these values we find 
1-2), 



core — ^- t *-corc 

(taking In 



^(coll) 
' Hydra 
^-j(Kolmo) 
' Hydra 



= (0.93 x 10~ 9 ) gl v. 
= (0.90-17.9 x 10 



2 

214, 
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510^214 ' 



If, on the other hand, we consider a specific example of 
a non-cooling core cluster such as the nearby Coma clus- 
ter of galaxies, we find little change to the chameleon- 
photon conversion rate. The IC magnetic field at the 
centre of the Coma cluster was determined from Fara- 
day RMs of an embedded source by Feretti et al. [32| . 
They found a uniform component of strength 0.2±0.1 fiG 
coherent on scales of 200 kpc and a tangled component 
of strength 8.5 ± 1.5 ^G coherent on much shorter scales 
of 1 kpc extending across the core region of the cluster 
(L co te ~ 198 kpc). The electron density at the centre of 
the cluster is 4 x 10 -3 cm -3 . We approximate L by L COIe , 
and find 



V, 



(cell) 
Coma 
(Kolmo) 
Coma 



(1.4 X 10- 9 ) ffl >2 2 14, 



0.64-12.7 x 10" 



2 5/3 
510^214 ■ 



(16) 
(17) 



We note that although the Coma cluster core has a 
weaker magnetic field compared to Hydra A, it also has 
a lower electron density which enhances the chameleon 
effect. Hence the two factors approximately balance out 
to give a similar magnitude of the CSZ effect for both 
Hydra A and Coma. 



B. The Coma Cluster 

The nearby Coma cluster of galaxies is unique in hav- 
ing detailed measurements of both its magnetic field 
structure and the SZ effect. Using SZ measurements 
along lines of sight through the core of the Coma clus- 
ter, we are able to constrain V-y^cj, for the Coma cluster. 
The SZ effect towards the centre of the cluster has been 
measured over a range of frequencies by various probes: 
OVRO, WMAP and MITO. These results have been as- 
similated and analysed in [33[. The frequency measure- 
ments of the change to the CMB background temperature 
are summarized in Table |H] and plotted in Fig O 

The predicted thermal SZ effect was given in Eq. (fl~3"l) 
and depends on the optical depth and the electron tem- 
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FIG. 2: SZ measurements of the Coma cluster from OVRO, 
WMAP and MITO. The solid line is the predicted ther- 
mal SZ effect assuming a maximum likelihood optical depth, 
ro = 4.7 x 10 -3 . The dashed line is the upper bound to 
the maximum likelihood fit for the combined thermal and 
chameleonic SZ effects. 



perature in the plasma. For the Coma cluster atmo- 
sphere, kT c ~ 8.2 keV. The optical depth t is gener- 
ally treated as an unknown parameter and inferred from 
the intensity measurements using a maximum likelihood 
analysis. Including the chameleonic contribution which 
is assumed to scale as u>~ a , the predicted form for AT/Tq 
in the Coma cluster is, 



AT? 



red 



1.6 x 10~ 2 (xcoth (0 -4)to (18) 

+ (e- - 1) X— 1 (^) a " PcomaM, 



where x — oj/fceTo; ao = a if 
power spectrum model and ao = 



-2 < a < -1 in the 
—2 in the cell model; 



and 7 ? coma(wo) = Py^<f>(L,coo) for lines 

of sight through the centre of the Coma cluster. We 
maximize the likelihood, L, over the parameter space of 
r and 'Pcoma(wo) where 



(AT obs - AT 
21ogi = ~ 1 1 



npred 



(TC^Coma)) 2 



Here AT° bs are the observed values of AT at a frequency 
vi = oji/2n, and are their standard errors. AT^ Icd are 
the values of AT as given by Eq. (p~5|) . We define t 
and 'Pcoma to be the best fit values of To and 'Pcoma- 
Confidence limits are then estimated by assuming that 



X 2 (T0,^Coma) = -2 log 



L(T ,Pcoma) 
. L(t , Vcoma) 



follows a Xi distribution. We find that when u 
0.844 meV (v = 204 GHz), the constraints on V Co 
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are almost independent of olq. For all — 2 < ceo < — 1, 
the best fit value and estimated standard error of r 
is virtually independent of uq. Fig. [3] shows the 68%, 
95% and 99.9% confidence limits on To and Pcoma^o = 
0.844 me V). Maximizing \ 2 with respect to "Pcoma; 
gives to = (4.7 ± 1.0) x 10~ 3 . Maximizing with re- 
spect to To gives the following 95% confidence limit on 
7> Coma (204GHz): 

^Coma(204GHz) < 6.2 x 10~ 5 (95%). 

The dashed line in Fig. [5] shows the best fit total 
SZ effect for the 95% confidence upper-bound when 
"Pcoma(204GHz) = 6.2 x 10~ 5 , with a = -2. It is clear 
that the preference towards smaller values of Pcoma is 
due mostly to the measurement points at 143 GHz and 
272 GHz which have the smallest error bars. A very sim- 
ilar picture is seen for other values of cvq. From Eqs. ([IT 



Confidence Limits from 
Coma Cluster SZ Measurements 




FIG. 3: Confidence limits from Coma cluster SZ measure- 
ments on the optical depth, to, and the photon to scalar con- 
version probability, Pcoma, at 204 GHz. 

and (fl"7|). we have 

P££(204 GHz) = (1.3 x 10-^) 5 ? (f^)' 

/ 4xlO~ 3 cm~ 3 N \ 2 / L \ 
\ n J \ 198L C oh J ' 



(204 GHz) = (0.59-11.8 x 10~ 7 ) 5 2 

f 4xl0- 3 cm- 3 > r 5/3 ( L \ ( L coh \ 
\ n e J \^ 198 kpc J \^lkp C y 



-2/3 



Fixing Brand = 8.5 fiG, L = 198 kpc, L con = 1 kpc and 
n c = 4 x 10 cm as suggested by observations, we 
find that the 95% confidence upper bound on T-'coma is 
equivalent to 

, 9 ^ cll) <2.2xl0- 8 GeV- 1 
in the cell model for the magnetic field, or 

5 (Koimo) < (7 2 _ 32 5) x iQ-io GeV" 1 , 



if turbulent Kolmogorov power spectrum is assumed. We 
note that the constraint g c g < 7.2 x 10~ 10 GeV , which 
corresponds to In = 1 + (<5n 2 /n 2 ) w 2 and the power 
spectrum being assumed to be of Kolmogorov form, is 
equivalent to M cS > 1.1 x 10 9 GeV. This was the con- 
straint on M e ff found by limiting the chameleonic produc- 
tion of starlight polarization in Ref. 0]. It is clear that 
whilst the constraint on V 1 ^$ at 204 GHz is reasonably 
independent of the magnetic field model, the equivalent 
constraint on g c g is not and depends strongly on the form 
and magnitude of the fluctuations in the magnetic field 
and the electron number density. 



C. SZ Profiles 

The above analysis for the Coma cluster involves mea- 
surements of the average intensity decrement towards the 
centre of the cluster at multiple frequencies. Higher fre- 
quency measurements are more successful in constrain- 
ing the conversion probability than low frequency ones, 
since Aicsz grows with frequency. Potentially better 
constraints could be found by also considering the radial 
profile of the intensity decrement. In §III| we discussed 
the scaling of the magnetic field and electron density with 
radius: 



1 + 



-3/3/2 



2\l/2 -n 

oc n', 



(B 2 ) 

where r c is the core radius and we have assumed a 0- 
model for the gas distribution. The standard thermal 
SZ effect is proportional to the integral of n e along the 
line of sight, while the chameleonic SZ effect depends 

1/2 

on n e , (B 2 ^ and £ C oh- As one moves away from the 
cluster core, simulations show that the coherence length 
increases, however the radial scaling L co h is uncertain. 
To calculate the radial dependence of the chameleonic 
SZ effect, we must specify how L con depends on r. We 
make the simple ansatz that, at least inside the virial 
radius, £ con oc n c 7 ^ 3 ' 3 for some 7, so that 



^cohO 



1 



7/2 



where £ C oho is the coherence length in the centre of the 
cluster. Based solely on dimensional grounds, estimated 

— 1/3 

values for 7 are: 7 = p so that L co h oc n c , or 7 = 1 so 
that L co h oc r for large r. The simulations of Ref. [24j . 
suggest that for large r, L con roughly scales as r inside 
the virial radius, r vlI . 

For a constant (<5B 2 ), n e and Z/ C ohj we found that, 

where ao = aif— 2 < a < — lin the power spectrum 
model, and cvq = —2 in the cell model. When (<5B 2 ), n c 
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FIG. 4: Chameleonic (red dot-dashed line for a — — 2 and solid blue line for a — —5/3) and thermal (dashed black line) SZ 
intensity decrement profiles for a simulated cluster at six different frequencies (taken to be the same as the frequency bands in 
the Coma cluster measurements). The simulated cluster has P 7< _ > ^ ) (204 GHz) = 3x 10~ 5 and To = 5 X 10~ s through the centre 
of the cluster, r c — 100 kpc and r v ir = 2 Mpc. We assume r\ — 0.9 and f3 = 7 = 1. 



and L C oh vary slowly along the line of sight we should 
change this to 

TV^cx J dzLlt^{z){5B 2 {z))nT{z)dz. 

We define r± to be the projected distance from the cluster 
centre for a given line of sight. We also assume that the 
turbulent cluster magnetic field only extends out to some 
characteristic radius r max , and then dies off very quickly 
due to i co h growing very quickly or (<5B 2 (z)) decreasing 
much faster than n2- We discuss reasonable values for 
r max below. With our assumed models for the radial 
dependence of the coherence length, magnetic field and 
electron density, we have 

V 1 ^,p{r ± ) oc L COIC ^1 + I x {r±;p,r max ,r c ), (19) 

where p = 1 + (1 + 00)7 — 3/3(2/7 + ao), L COIC — 2r c and 

roc (p-1) 

I x = dx (1 + x 2 ) 2 w{x;x mliK ). (20) 
Jo 

We have introduced a window function, w(x), to impose 
the rapid decay of magnetic turbulence for r > r max . A 
simple model for w(x; x max ) would be that w = 1 for 
r 2 < r 2 lax and otherwise. This corresponds to w = 
H(x ma x — x) where H is the heaviside function, and 



^max 

r 2 + r\ 



Alternatively, we could impose a smoother decay of the 
turbulent magnetic field by taking w — exp(— r 2 /r 2 lax ), 
which corresponds to 

W (,)=expf- Ii±4±^ ). (21) 

In what follows we approximate w(x) by this latter form. 

All that remains is to determine a reasonable estimate 
for the largest radius to which the turbulent cluster mag- 
netic field extends, r max . We note that whenever p > 0, 
the dominant contribution to the integral in Eq. (|19[) will 
come from the largest values of x, i.e. from r ~ r max . 
The value of r max will, therefore, play an important role 
in such scenarios. Assuming [3 — 7 = 1, in the cell mag- 
netic field model (ao = —2) we have p cc ii = 6(1 — rj) and 
so p ce ii > for rj < 1 , while if we assume a Kolmogorov 
spectrum with ao = —5/3 we have PKni mn = 6(8/9 — rj) 
and so p cc \\ > for 77 < 8/9. In Ref. [3J], the quan- 
tity n7 rms t age is introduced where i agc is the age of the 
Universe and tn rms is the root mean squared value of the 
vorticity' m — V x v, where v is the peculiar velocity. 
In Ref. [34j | it is estimated that when w Ims t ase is greater 
than O(l), turbulence has had time to develop. In a viri- 
alized structure we have the following order of magnitude 
estimates: m ~ 0(v(r)/r) andf 2 (r) ~ 0(M(r)/r) where 
M(r) is the mass inside a radius r. Thus, 

where p(r) is the average density inside a radius r and 
p c is the cosmological critical density today. Beyond the 
virial radius, we expect that vd 2 < 0(M(r)/r 2 ), as the 
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FIG. 5: The effective path length L through the magnetic region. The leftmost plot shows p values corresponding to r\ — 0.9 
and rj = 1 in the Kolmogorov (a = —5/3) and cell (a = —2) magnetic field models when = 7 = 1. These are the values of 
r\ suggested by observations and simulations. The rightmost plot shows p values corresponding to 77 = 1/2 and 77 = 2/3 again 
in the Kolmogorov and cell magnetic field models with /3 = 7 = 1; these are the values of rj suggested theoretically. In the 
Kolmogorov model r\ — {1/2, 2/3, 0.9, 1} corresponds to p — {7/3, 4/3, —2/30, —2/3}. In the cell model the same values of 77 
give p= {3, 2, 6/10,0}. 



vorticity acts to slow down the collapse. Roughly then we 
expect w rms t agc to be smaller than a few when r 3> r vir , 
and larger at smaller radii. We choose r max so that the 
chameleonic SZ effect dies off quickly for r > r v „. With 
the choice Gaussian window function given in Eq. (|21[). 
a suitable choice is r max w r v ir- 

We have the following approximate scalings for the 
thermal and chameleonic SZ effects: 



AIsz 
To 

AXcsz 



1 + 



1 + T 4(rx;p,rvir,r c ), 



where T is given by Eq. (f^U]) and we take the Gaussian 
window function given in Eq. (|2ip to impose the decay 
of the turbulent magnetic field for r > r max r v ; r . For 
rj_ *C r V i r , the chameleonic contribution to the SZ effect 
will decrease less rapidly than the thermal one when, 

p + 3/3 - 1 = (1 + a ) 7 + 3(1 - ao - 2rf)P > 0. 

For 7 ps f3, and 77 < 1 this holds for all —2 < ao < — 1- 
Hence, even if the thermal SZ effect dominates near 
the centre of the cluster, measurements of the SZ ef- 
fect at radii r c <C r <C r v - 1T might be able to detect the 
chameleonic effect. 

In Fig. [4] we plot the predicted profiles for the 
chameleonic and thermal SZ effects at different frequen- 
cies for a model cluster, with 77 = 0.9 and (3 = 7 = 1. We 
assume that the cluster properties are similar to those 
of the Coma cluster and assume V-y^,^ = 3 x 10~ 5 at 
v = 204 GHz for a line of sight through the centre of 
the cluster. This is roughly the la confidence limit de- 
rived earlier from Coma cluster observations. We take 



r c s» 100 kpc, and r max s» r v ; r w 2 Mpc. The opti- 
cal depth for a line of sight through the centre of clus- 
ter is taken to be To = 5 x 10~ 3 . The black dashed 
line shows the usual thermal SZ effect, the dot-dashed 
red line is the chameleonic SZ effect for ao = —2 (cell 
model) and the solid blue line is the chameleonic SZ ef- 
fect for ao = —5/3 (Kolmogorov turbulence model). In 
general, the smaller the value of 77, the more slowly the 
chameleonic SZ effect decays with distance from the clus- 
ter centre. We note that for 30 GHz < v < 300 GHz, a 
chameleonic SZ effect of this magnitude dominates over 
the thermal SZ effect for r± > 500 kpc. At frequencies 
around 150 GHz < v < 275 GHz, the chameleonic SZ ef- 
fect is dominant for r± > r c . Measurements of the CMB 
intensity decrement associated with clusters over such a 
frequency range and for r > r c could be particularly use- 
ful in constraining the properties of chamcleon-like theo- 
ries. 

There has been some suggestion of a discrepancy be- 
tween X-ray emission data as measured by ROSAT, 
which determines the electron density distribution, and 
the SZ intensity decrement as measured by WMAP [35[. 
In the analysis the expected signal derived from X-ray 
measurements is fitted to the observed SZ effect at large 
radii. They conclude that the observed intensity decre- 
ment at the centre of the cluster is then less than the 
predicted SZ effect. However a possible alternative is the 
existence of a chameleon-like scalar field, which would 
cause an enhancement in the SZ effect at large radii. The 
WMAP data includes measurements at 94 GHz for which 
the chameleon effect could be significant. 

In W Al we estimated "P 7 <_>0 in terms of the properties 
of the cluster magnetic fields, and approximated the path 
length, L, through the magnetic region for a light ray 
passing through the cluster as being equal to twice the 



core radius, i.e. L 



2r c . However, the relatively 
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flat radial profile (for r < r v i r ) of the chameleonic SZ 
effect derived above suggests that L is generally much 
larger than L coro , and for p > is better approximated 
as being 0((r vil /r c ) p L COIC ). Using the model introduced 
above for the radial behaviour of the electron density, the 
magnetic field and L co h, we find that the effective path 
length is given by: 



L — L COIC I x (0; p, r n 



cell and power spectrum models we find, 



(Qq) 

<Qu) 
(Qv) 



= (Q-Qo) = Qocos20 r 
= (U-Uq) = Q sin26l rc 
= (V-V ) = 0, 



and B r y c s = S r °s sin rcg define 



where Bf*> = B lc s cos 

9 TCg . Detailed derivations of these formulae are presented 
in Appendix A. In the cell model, 



It follows from Eq. $2$ and Eq. (51]) that L/L corc 
depends on two parameters p and r V i r /r c . In Fig. 
we plot L/L COIC against r v - n -/r c for different values of 
p. Reasonable estimates for r\ are rj = 1/2, 2/3, 0.9 
and 1. 77 = 1/2 and rj — 2/3 are based upon theo- 
retical expectations, the former being for equipartition 
of magnetic and thermal energy and the latter being if 
the magnetic field is frozen in. Observations and simu- 
lations suggest that these values are unrealistic and in- 
stead indicate rj = 0.9 and rj = 1 respectively. With 
(3 = 7 = 1 and rj = {0.9,1}, p = {-2/30,-2/3} in the 
Kolmogorov model, and p — {6/10,0} in the cell model. 
We plot these values of p for 5 < r v i r /r c < 40 in the 
leftmost plot of Fig. [5] We can see that in all four cases 
icorc < L < llL C ore for the range of r V i r /r c considered. 
Since P 7 „^, oc g^ s , this corresponds to no more than a 
factor 3 or so enhancement in sensitivity to <? c ff. If in- 
stead we consider rj = {1/2,2/3} then p = {7/3,4/3} in 
the Kolmogorov model and p = {3, 2} in the cell model. 
We plot these values of p in the rightmost plot of Fig. [5] 
In these cases, we see that L 3> L cmc , and as expected 
L/L cole ~ 0((r vil /r c ) p ). Since rj = 1/2 and rj = 2/3 do 
not appear to be supported by observations, and because 
the structure of cluster magnetic fields at radii r ^> r corc 
is far from certain, we do not view values of L > £ C orc 
as being particularly realistic. The more realistic scenar- 
ios are those shown in the leftmost plot in Fig. \5\ where 
1 < L/L cmc < 11. 



VI. MODIFICATIONS TO THE CMB 
POLARIZATION 

In addition to modifying the intensity of light that has 
passed through magnetic regions, we found in SjlT] that 
the presence of a chameleon-like scalar field also alters 
the polarization of the light. We define Qo, Uq and Vq to 
be the fractional intrinsic CMB polarizations. We found 
that after mixing with the chameleon in a magnetic re- 
gion, 

Q = Qo + Qq, 
U = Uo + Qu, 
V = V + Q v . 



1 / 2B rC giu 



Qo — 77 l T7 r 

2 \M cS m z cS 



sin 



2 A, 



whereas in the power spectrum model, 



Qo 



1 / 2B T , 



4 \M cff 77l^ 

1 / 2B rC gLu 



r 3 



4 \M oS m 2 cS 



cos 



Vff 



L 



2oj 



In both cases Qo = P™ 8 where p re & is the contribution of 
the regular magnetic field to the conversion probability. 
In addition to the contribution from mixing in the regu- 
lar magnetic field along a given line of sight, the Stokes 
parameters will also pick up a contribution from the tur- 
bulent magnetic field. The magnitude and direction of 
this contribution is determined by the direction and size 
of the magnetic field fluctuations along the line of sight, 
and so vanishes when averaged over many lines of sight. 
Along any given line of sight it would be observed as an 
essentially random contribution, and would be present in 
the variance of the Qi, given by of (A, A): 



o 2 (A l5 A 2 



(e*(Ai)Gi(Aa)) - (Qi(Ai)) (Qi(A 2 )) 



A detailed calculation of these <Ji is given in Appendix 
A. We find similar results for both the cell and power 
spectrum models. In the latter we have, 



°"q(A, A) 
°u(A, A) 
^v(A,A) 



2 
1 

2 

2 ^ 



-|2 . 

2" 

2' 
1 



7^2ocos40 re 



Hence, looking along a single line of sight at a single fre- 
quency, we would expect \Q\, \U\, \V\ - 0(P 7 ^/\/2), 
since 'P 7 <->^ is dominated by the contribution from the 
random magnetic field (V^,/, 3> V 1CS — Qo)- The best 
upper bound on V-y^tj, from the Coma cluster SZ mea- 
surements was, 

P 7 ^(204GHz) < 6.2 x 10~ 5 , 



We use angled brackets to denote an average over many 
lines of sight at fixed frequency, i.e. fixed A. In both the 



along lines of sight passing though the Coma cluster core. 
Along such lines of sight, we might expect chameleonic 
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contributions to Q, U and V that are as large as a 
few x 10~ 5 . If this were the case then the magnitude of all 
three Stokes parameters would be roughly 10 — 100 fiK, 
which is orders of magnitude larger than the measured 
0.1 — 1/iK level of Q and U and the even lower bound 
on the circular polarization, V. One might imagine, 
therefore, that CMB polarization measurements could 
greatly improve the upper limit on V 7 ^^>, say placing 
it around the 10~ 7 level. We found in fV Bl that in the 
Kolmogorov model for the power spectrum of magnetic 

fluctuations, Vco^ could be almost as large as 10~ 6 
with gio ~ 0(1). An upper limit of 10~ 7 could there- 
fore constrain M e g > 10 10 GeV. However, there are two 
important caveats of which we must take account and 
which limit the ability of realistic experiments to detect 
the chameleonic contribution to the CMB polarization. 

Firstly, the random effect (whose variance is given by 
the of) vanishes when one averages over many lines of 
sight which have taken different paths through the ran- 
dom magnetic regions. The coherence length, L con , pro- 
vides the typical length scale of the random magnetic 
regions. We therefore expect that roughly parallel light 
rays with a perpendicular separation d 3> L co h will have 
taken different paths through the magnetic regions. If 
the cluster in question is a distance, D c lust 3> d, from us 
then light rays with a separation, d < L co h, will subtend 
an angle < (9 co h when observed on the sky, where 

. 10800 L coh 
fcoh = — aremm. 



D, 



lust 



As an example we consider the values for these quanti- 
ties for the Coma and Hydra A clusters: (D c iust, ^coh) ~ 
(100 Mpc, 1 kpc) and ~ (50Mpc, 4kpc) respectively. 
This results in # co h ~ 0.03 arcmin for Coma and # co h ~ 
0.27arcmin for Hydra A. In order to resolve the random 
chameleonic contribution to the CMB polarization, a sur- 
vey must have an angular resolution better than con . 
The polarization of the CMB in the direction of galaxy 
clusters has not so far been measured, but there are plans 
to incorporate a polarization detector at the South Pole 
Telescope, SPT-pol, which is expected to be completed 
in the next few years (36|. SPT is searching for clusters 
on the sky by exploiting the SZ effect. Each pixel of de- 
tectors at SPT represents one arc-minute, and the aim is 
to have 10 /iK sensitivity per pixel. For a cluster such as 
Coma where L C oh is relatively small, this one arc-minute 
angular resolution is not enough to resolve the random 
chameleonic contribution to the polarization. Even for a 
cluster such as Hydra A where L co iJ D c \ ust is larger, it is 
still a factor of four larger than # co h. In general, resolving 
the random chameleonic contribution to the Stokes pa- 
rameters will require an ability to resolve the polarization 
signal on scales smaller than a few arc-seconds. 

Even if one has the required angular resolution, there 
is still an additional technical limit on the ability to de- 
tect the random chameleonic modification to the Stokes 
parameters. The expressions given for the of assume 
that measurements are made at precisely one frequency. 



Any realistic measurement will, however, not be entirely 
monochromatic but represent an average over some fre- 
quency band. Suppose one makes a measurement at 
a nominal wavelength, A, which is actually an average 
(with some window function) over a range A — A e 
(—SX,8X). This corresponds to the value of A running 
between A — 6 A and A + SA where 6 A = A(SX/X). The 
measured variance in the Stokes parameters will then be 
given by, 



1 



(2SA) 2 



A+SA i-A+SA 

dAx / dA 2 ^ 2 (Ai,A 2 ) 

ASA JA-SA 



A full calculation of the of (Ai, A 2 ) shows that the domi- 
nant term in it drops off as l/(Ai — A 2 ) when |Ai— A 2 | ^> 
1. It follows that when <5A 1, of is suppressed by a 
factor 0(1/ 8 A 2 ) relative to its value when 8 A = 0. Now, 



8A 



(4.3 x 10 3 ) ( 



io- 



L 



8X 



100 kpc / \1 /im 



Typically for CMB radiation A ~ 1 mm-1 cm, and so 
to achieve 8 A < 1, we would require 8X/X < 10~ 6 — 
10~ 7 which is many orders of magnitude smaller than the 
typical spectral resolution of CMB experiments, 8X/X ~ 

io- 4 -io~ 2 . 

It is clear then that both the angular and spectral res- 
olution of detectors would have to improve significantly 
before the random chameleonic contribution to the CMB 
polarization could be detected. The smaller regular com- 
ponent, however, does not suffer from such issues, al- 
though unfortunately it does not source the (more un- 
usual) circular polarization in the CMB. In the cell model 
the contribution from the regular magnetic field is a fac- 
tor of 1/N smaller than Pj^. In the power spectrum 
model, however, electron density fluctuations result in 
some enhancement of the regular magnetic field conver- 
sion probability. The regular contributions to both Q 
and U are generally dominated by the term in Qo com- 
ing from the electron density fluctuations: 



Qo 



B 10g L 

16A/> C 2 



Assuming that fc 2 P N = 2(2n) 1 / 3 C^k a for k > L"^, we 
found in mil AI that for a Kolmogorov spectrum with 
a = —5/3, 



2(2tt) 1 / 3 C 2 i 



(0.27-0.45)n 2 (^ 



I'll— 2/3 



Hence in the Kolmogorov model, 



-2/3 



n 1(1 -2^UN - 1)-B IC g LL coh ,-5/3 

IVI eS 

w (0-1.5 x lO- 9 )g 2 w 4[lB 2 llJtG L 200 

-5/3 / .coh N-2/3 
"■0.01 l^lkpc/ ) 

where as before 7 N = 1 2, g w = 10 10 GeV/M cf j, etc. and 
we define Bi^q = _B rog /l/iG. For the Coma cluster with 
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L w 198 kpc, n c w 4 x 10" 3 cm- 3 , £ rog w 0.2 /iG and 
icoh ~ 1 kpc this gives, 

2„ Coma ~ (0-2.9) x 10- w g 2 10 4(l 

whereas for Hydra A with L w 260 kpc, n w 10~ 2 cm -3 , 
L C oh ~ 4 kpc and S rcg w 6 /zG we have, 

Q« ydra ~ (0-3.0) x lO-Vo^t 

For the chameleon-induced linear polarization to be at a 
similar level to the intrinsic CMB polarization one would 
require Qq ~ fewx 10 -7 . For the Coma cluster this would 
require g w > 27 corresponding to Pcoma(214 GHz) ~ 
few x 10~ 4 , which is strongly ruled out by our analysis 
of the SZ effect. However in Hydra A, where the large 
scale magnetic field is much stronger, one could achieve 
<2o ~ few x 10 -7 with only g%o sa 2.6 which is a factor 
of 4 smaller than the lowest upper bound on gio found 
previously from the Coma SZ measurements. It is feasi- 
ble, therefore, that future searches for CMB polarization 
in the direction of galaxy clusters could, in some cases, 
reveal an enhancement of polarization at higher frequen- 
cies. However this would only be expected to occur for 
clusters, such as Hydra A, with suitably large regular 
magnetic fields, O(10/iG). 



VII. CONCLUSIONS 

The chameleon was first introduced to explain the pos- 
sible scenario of a very light scalar field existing in the 
cosmos with a gravitational strength coupling to mat- 
ter and yet not being detected on Earth as a fifth-force. 
The theory has since been extended to allow supergrav- 
itational strength couplings and different forms of the 
matter coupling. Recently, there has been interest in pre- 
dicting the astrophysical implications of a direct coupling 
between chameleon-like fields and the photon [3, y, 0]- 
This analysis applies not only to the chameleon model 
but also other chameleon-like theories such as the Olive 
Pospelov model designed to explain a variation in the 
fine-structure constant, a em - 

In this paper we have specifically considered the ef- 
fect of a chameleon-like scalar field on the CMB intensity 
and polarization in the direction of galaxy clusters. The 
chameleon mixes with CMB photons in the magnetic field 
of the intracluster medium. The resulting modification 
to the CMB intensity could be observed as a correction 
to the existing Sunyaev-Zel'dovich effect. Measurements 
of the CMB polarization in the direction of clusters are 
still waiting to be made, but the chameleon effect could 
potentially be significant in the polarization foregrounds 
as well. 

Predictions of the chameleon effect are strongly depen- 
dent on the assumed structure of the intracluster mag- 
netic field. In this work we have considered two differ- 
ent models for the magnetic structure: a simplistic cell 
model in which the cluster properties are uniform but 



the magnetic field has a random orientation in each of 
the N cells along the light path; and a potentially more 
realisistic power spectrum model in which fluctuations in 
the magnetic field and electron density are allowed over 
a range of length scales, with the small-scale power spec- 
tra taken to correspond to three-dimensional Kolmogorov 
turbulence. 

The conversion probability, 'P 7 ^0, in the cell model for 
magnetic fluctuations is given by, 

V^(L^) = (3.2 xlO- 10 )^^^ 

i20o(il°kpc) ^214) 

while in the Kolmogorov turbulence model it is predicted 
to be slightly larger with 

5/3 



7^ „2 n 2 

lO^G^O.Ol 



p( Kolmo )(L,z/) = (2.4-3.8 x 1CT 7 ) g\ B 

5/3. fj- coh y 2/3 / /5(2J N -1) \ 
^214^200 ^1 kpc J y 12 )■ 

The resulting change to the CMB intensity contributing 
to the SZ effect is given by, 



AL 



csz 



Io 



(L). 



The strong dependence of the conversion probability 
on frequency allows us to compare measurements of the 
SZ effect across a range of frequencies, with our predic- 
tions for the combined thermal SZ and chameleonic SZ 
effects. A maximum likelihood analysis of the SZ obser- 
vations of the nearby Coma cluster leads to a constraint 
on the chameleon-photon conversion probability: 

^Coma(204GHz) < 6.2 x 1CT 5 (95%CL). 

This can be converted into a constraint on the chameleon- 
photon coupling strength. However this is strongly de- 
pendent on the magnetic fluctuation model assumed and 
also on the properties of the Coma cluster. Assuming the 
current available data on the Coma cluster properties [HJ 
we find, 

5 ^ cll) <2.2xlO- 8 GeV- 1 
using the cell model, and 

5 (Kolmo) < (7 2 32 5) x 10 -10 GeV -l ; 

assuming a Kolmogorov turbulence model. We note that 
other clusters such as Hydra A, a good example of a 
cooling-core cluster, have approximately similar conver- 
sion probabilities, since although the observed magnetic 
field strength at the centre of the cluster is much greater 
(enhances chameleon effect) the electron density is also 
much higher (supresses chameleon effect). 

The model dependence of the effective coupling 
strength leads us to reevaluate the constraints given 
in Ref. Q. Assuming a cell model for the magnetic 
field fluctuations the constraint on the coupling strength 
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from starlight polarization was found to be g c g < 9.1 x 
10~ 10 GeV . From our analysis we have found that the 

conversion probability increases by a factor of (A/AT) 1 
if a Kolmogorov turbulence model is assumed instead. 
This increases the sensitivity to g e s by approximately 

(A/N) 1/6 . In our galaxy, A/N ~ 50-100, for the ob- 
jects considered in Ref. [J]. Thus we expect to increase 
the constraint to approximately g e g < 4 x 10~ 10 GeV _1 
if a Kolmogorov turbulence model is used in the analysis 
of the starlight polarization data. We note that the cur- 
rent constraints coming from starlight polarization are 
tighter than those from SZ measurements. However for 
the starlight polarization the data is predominantly lim- 
ited by intrinsic scatter and improvements are likely to 
be statistical, while there is the possibility of SZ observa- 
tions away from the cluster core leading to even better SZ 
constraints (see below). It should be mentioned, however, 
that there is significant uncertainty in the SZ constraints 
coming from a lack of detailed knowledge about cluster 
magnetic fields on scales ~ fc~; t , whereas the starlight po- 
larization analysis rests on knowledge of the power spec- 
trum of the galaxy magnetic field and electron density 
fluctuations which are much better known. 

In addition to the frequency dependence of the 
chameleon effect, it is possible to exploit the depen- 
dence on magnetic field strength and electron density to 
examine the radial dependence within the galaxy clus- 
ter. The magnetic field and electron density are ex- 
pected to decrease with radius. We have assumed a 
simple scaling for the electron density following a (3- 

( 2 \ ~3/3/2 

profile, n e — hq I 1 + ^ ) , with B oc n^. In Fig. 

5] we plotted both the thermal SZ radial profile and 
the predicted chameleonic radial profile for a hypothet- 
ical cluster with chameleon-photon conversion probabil- 
ity "P 7 <^(204 GHz) = 3x 10~ 5 through the centre of the 
cluster. We found that the chameleonic effect generally 
dominates towards the edges of the cluster. The sugges- 
tion of a discrepancy between the inferred SZ effect from 
X-ray emission data and the WMAP SZ measurements 
[HI could potentially be attributed to a chameleon, and 
this requires further investigation. In any case, future 
measurements of the SZ profile at higher frequencies 
would lead to even better constraints on the chameleon 
model. 

Finally we looked at the modification to the CMB po- 
larization. We found that the contribution to the modi- 
fication from chameleon-photon mixing arising from the 
random component of the intra-cluster magnetic field, i.e. 
that tangled on scales much smaller than the size of the 
cluster, is severly limited by the angular and spectral res- 
olution of the detector. We would require angular resolu- 
tion on the scale of the magnetic domains in the cluster, 
approximately one-hundredth the size of the cluster, and 
frequency resolution of the order 8X/X < 10 -6 before 
any effect could be detected in the CMB polarization. 
The contribution arising from the regular component of 
the magnetic field is more significant and potentially de- 



tectable. So far polarization measurements of the CMB 
in the direction of galaxy clusters have not been made. 
The upcoming SPT-pol experiment [3(| will change this. 
Assuming a Kolomogorov turbulence model, we predict 
the change to the fractional CMB linear polarization: 

AQ ps AU « (0-1.5 x lO- 9 )g 2 w 4(lB^ G L 200 

-5/3 (T coh N-2/3 
''0.01 l^lkpcj i 

which for some clusters is approaching the level of the 
intrinsic CMB polarisation. 

We conclude by mentioning two points which at first 
may appear to be potential sources of error in our pre- 
dictions. Firstly we remind the reader that the magnetic 
field strengths of most galaxy clusters are determined 
from Faraday rotation measures, looking at the degree to 
which the polarization plane of some background source 
is rotated due to the presence of the cluster magnetic 
field. It is feasible that the presence of a chameleon 
could influence the change to the polarisation angle and 
hence the inferred magnetic field strength of these clus- 
ters. However the Faraday RMs generally come from 
low frequency measurements around 10 GHz for which 
the chameleon interaction is greatly suppressed. There- 
fore we believe it is realistic to ignore the contamination 
when considering the cluster magnetic field strength. 

Secondly, in these predictions we have neglected the 
contribution from our own galactic magnetic field. In Q 
the bound on the chameleon-photon conversion probabil- 
ity in our own galaxy coming from starlight polarization 
was found to be < few x 10~ 6 at u> = 1 eV. As- 

suming the probability scales as lu 5 ^ 3 , then the contribu- 
tion at CMB frequencies from the galactic magnetic field 
would be less than O(10 -10 ), which is negligble. 
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APPENDIX A: PHOTON-SCALAR MIXING 

In the weak mixing limit, we saw in JET] that the 
photon-scalar mixing can be described by a function 
Ai(z) where i — {x, y}, and 

Mz)= / 2 d S |Me^W. (Al) 

JO ^eff 

We assume that B has both a regular component, B rog , 
which we take to be constant, and a random component, 
(5B . We also found in SjTT] that the probability of a pho- 
ton converting into a scalar field, and the induced 
contributions, Q qi Q u and Q v to the Stokes parameters 
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Q, U and V, are given in terms of the A t thus: 
?W*) = \ {\A x {zt + \A y {zt), 

Qq(*) - i(|^(z)| 2 -|^(z)| 2 ), 

Q u (z) = Re(A* x (z)A y (z)), 
Q v (z) = lm(A* x (z)A y (z)). 

We now evaluate these quantities in the cell and power 
spectrum models for the magnetic field fluctuations when 
z = L, where L is the path length through the magnetic 
region. 



1. The Cell Model 

In the cell model for the magnetic field we assume that 
along the light path, < z < L, the components of the 
random magnetic field <5B perpendicular to the line of 
sight are given by 

(SB) X = B rand COs6>„, 

(SB) y = Brand sin 6»„, 

in the region (n — l)L C oh < z < n_L con , where L con is 
the coherence length and 6 n ~ U(0,2n]. We define N = 
L/L co h with n running from 1 to N. We take -B ran d to 
be fixed and assume a constant value of n c = n c . Thus 
A(x) — Ax/L, where 



and m 2 ff ~ — = — (47ra om n e /m ) 2 in the cluster at- 
mosphere. It is now straightforward to evaluate Ai and 
we find, 



ML) 
A rC g 

5A 



(A r0 g)i 



M eff m 2 ff 



(SA) U 

A 

sin A, 



M cS m l cS 



9 sm 



A 



where n v ' v = 

BrogCOS^rog and (B re g)y 



(«) = 



(cos #„,sin 9 n ,0) 

B r( ,„ sin 9 rc „ i wc find 



JV-l 

53 n (r) e^ 

r=0 



Taking (B r0 g) a 



rog 

v 7 ^4,(L) = r TCS + r cmss + -p rand , 

Qi(L) = Q? s + Qf oss + Qr d , 
where the regular components are given by: 

1 

2 



■prog 



Q rcg 



o reg 



cos 28 IC 



sin 2 A, 



2 

sin 29 r 



2B rcg u) V 

2B rcg uj 
M eff m 2 ff 



sin 2 A, 



sin 2 A, 



0. 



across 



Similarly the cross-terms are given by: 



2^ V f A 

— — BregBrandSinAsm — 

M oS m* s J 6 

JV-l s 

53 cos(# r - 9 Icg ) cos ( — (2r + 1 - JV) J , 
2w \ 2 /A 

— 2~ -Breg Brand Sin A Sin — 

JV_1 /A \ 

53 cos (^ + ^og) cos I — (2r + 1 - N) J , 

r=0 ^ ' 

2w \ 2 /A 

— — -BregBrand sin A SHI — 

M eS m* s J 6 
53 sin(0 r + 6» rC g) cos ( — (2r + 1 - JV) J , 
2u \ 2 -r f A 

— 5- BregBrandSinAsm — 

MesmiffJ \N 

N -i /a \ 

53 sin(0 r - 6 rcg ) sin I -(2r + 1 - N) ) . 



Qcross 



Finally the terms due entirely to the random component 
of the magnetic field are: 



■>rand 



1 / 2B ra nd^ y . 2 fA 

2 ^M eff m 2 J Sm U 



JV-l JV-l 



53 53 cos(e r - e s ) cos 



s=() r=0 



Q; 



rand 



1 { 2B land u> y . 2 (A 

2 \M cS mlJ Sm U 



JV-l JV-l 



53 53 c ° s (^'- + 0«) c ° s 



=0 r=0 



2A 



2A 
"/V 



(r-s)\, 



Qrand 



1 / 2B ran( ±UJ \ . 2 ( A 

sm (iV 



2 VMcfrm 2 



cir 



JV-l JV-l 



53 53 sin(0 r + e 8 ) cos 



s=0 r=0 



2A 

jV v 



- s. 



rand 



1 / 2B ran dCJ \ . 2 / A 



2 lM cff m 2 



off 



53 53 sin (° r - 6s ) sin ( - s ) ) • 



s=0 r=0 



Taking an average over many lines of sight through the 
cluster is equivalent to averaging over the 6 n . Doing this 
it is straightforward to check that the average values of 



21 



■P 7 ^>0 and the Qi are given by, 



<*W = ' f^^) -rA 



2 V^effW 2 



cir 



N ( 2B rand w 



2 V^effTO^ ff 



A 

sm" [ — 

N 



(Qi 



We can also calculate the variance of the Qi over many 
lines of sight. Assuming |A|/JV ^> 1, N ^> 1 and remem- 
bering Brand ~ 0(B Icg ), we find 



i-p2 
2 ^ 



This assumes that measurements of the Stokes parame- 
ters, and hence the Qi, are exactly monochromatic. The 
observable erf can be greatly suppressed relative to the 
above estimate when, as is always the case in practice, 
measurements actually comprise an average of the ob- 
servable quantity over a frequency band about the de- 
sired measurement frequency. We discuss this further in 
~3l below. 



2. The Power Spectrum Model 

In the power spectrum model, we assume that fluc- 
tuations in B occur on all spatial scales, and we de- 
scribe the magnitude of the fluctuations on different 
scales by a power spectrum, Pb(&)- We also allow fluc- 
tuations in the electron number density, n c , which are 
described by another power spectrum, P^(k). We define 
n c (z) — n e (l + <5 n (z)), where n e is the average electron 
density along the path through the magnetic region. We 
also define, 



1 



S n (x)dx . 



Note 5 n (L)_= 0. With A, defined by Eq. JATJ, we define 
Z = (1 + 5 n (z))z and A = m 2 s L/4uj oc -UcL/Auj. It 
follows that A(z) = AZ/L. We assume that \S n (z)\ < 1 
along the majority of the path, and so Bi(z) w Bi(Z). 
We then have, 

1(X1 f L , BAZz) 

Ai(A) w / dZ 7 l±—±—— e —. 

To calculate the expected values of "P 7 ^0 and the Qi, we 
must calculate the correlation at different frequencies: 

A lJ (A 1 ,A 2 ) = (A l (A 1 )A*(A 2 )) 

where the correlation function is defined, 
BiiyfyBjdy + x)z) 



[l + 6 n (yz)(l + 5 n ((y + x)z) 



We split the magnetic field into a regular and turbulent 
component: B = B rog + 5B. Assuming isotropy of the 
fluctuations in B and n c , we rewrite thus, 

R\f(xz) = B^B^ 

■({l + 5 n (yz))- 1 {l + 5 n {{y + x)z))- 1 ) 



SB(yz) ■ SB((y + x)z) 



3 \(l + 6 n (yz)(l + S n ((y + x)z)) 



We note that isotropy of the magnetic fluctuations im- 
plies that the off-diagonal components of are sourced 
only by the regular magnetic field. We define a combined 
power spectrum P* ot (/c) by, 



/?;•/' f.rz) = I d 3 kP^ ot (k)e 27rlxk i . 



We then have 
L 2 



A 

where 



" " AM 2 S 



OO 



e <(Ai-A a ) J fc 2 P* ot (fc)^(^fcL;A 1 ,A 2 )dfc, 



K(irkL;Ax,A 2 



x , sin(27rfeis) sin(2A_(l - s)) 

ds = 

2irkLs A_ 

■ cos(2A + s), 



and we have defined 2A± = A2 ± Ai. This integral may 
be evaluated explicitly in terms of the functions: 

f x sin^/2) 
Gm{x) = / — — — dy, 



Si(x) 



y/2 

sin y 



V 



dy. 



We define 6Cin(x, y) — Cin(x + y) — Cin(a; — y) and 
<5Si(a;, y) = Si(x + y) — Si(x — y). We then have, 

cos2A_ 

4xK(x;A u A 2 ) = — [ < JCin(2Ai,2a;) 



2A_ 
-<5Cin(2A 2 ,22;)] 
sin2A_ 
2A_ 
+<5Si(2A 2 ,2x)] 



+ 



[5Si(2Ai,2x) 



We are concerned with CMB radiation, for which ui ~ 
10 _5 -10 -3 eV, propagating over distances of the order 
of 100 kpc through clusters. It is easy to see that for 
such a scenario |A| 3> 1. We evaluate K(x; Ax, A 2 ) 
in the asymptotic limit where |Ai|,|A2| 3> 1. When 
min(|Ai|, |A 2 |) > 1 and x < min(|Ai|, |A 2 |) we have, 



4xK(x;A 1 ,A 2 ) 



2x cos 2A_ 
AiA 2 

sin 2x cos 2A + 
AiA 2 : 
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and when x » min(|Ai|, |A 2 |), we have 

27rsin2A_ 2cos2A- 



4xK(x; Ai, A 2 ) 



2A_ 



We define fcj = \Ai\/wL. We assume that the fluctua- 
tions in B and n c are such that Pf° t {k) drops off faster 
than fc~ 3 for all k > /c* where fc* <C min(fci). This 
is equivalent to assuming that the dominant contribu- 
tion to R^fo) comes from spatial scales than are much 
larger than fcr . We expect that this dominant contribu- 
tion will come from scales of the order of the coherence 
lengths of the magnetic field and electron density fluctu- 
ations (Lb and Ln respectively) and so we are assuming 
that k^ 1 -C Lb, Ln- With this assumption, in the limit 
of large |Aj|, we have to leading order, 



L 2 cos2A_ 
2A!A 2 jo 
L cos 2 A 



k 2 P^°\k)dk 



+- 



4ttA 
L sin 2 A 



8A 



OA f 00 

--± / ksm(2-KkL)P£Hk)dk 
iA 2 J Q 

kPfi* (k)dk 



fci 



+ 



^ fc^ ot (fc)dfc 



Wc define, 



/>OG 

Wi J -(fc)= / ?P* ot (?)d g , 
Jk 



and we can then rewrite Aij as, 



. 2i A_ , ^ 4 Wl o; 2 . jR tot (0)cos2A _ 



eff"'eff 

^-«!f(^)co S 2A + 



'eff'eff 

Lsin2A_ 
~32MpT 



[Wii(fei) + W«(fe2)]. 



For simplicity, we assume that fluctuations in the mag- 
netic field and electron density are uncorrelated. Thus 



•((l + 5 n (y))- 1 (l + ( 5 n (y + x))- 1 ). 



This implies that, 

i?-(o)«(sr g sf s +^(' 5B2 ))(S)- 



We assume that fluctuations in n c follow an approxi- 
mately log-normal distribution, and hence 



n e\ _ / n c\ 3 - 7-3 



We assume, as is the case in the situations we consider, 
that L is much larger than the correlation length of the 
magnetic and electron density fluctuations. This implies 
that the only contribution to ^-^(Lz) which may result 
in a leading order contribution to Aij is, 

R*?(Li)~Br s Bp. 

Finally, we must evaluate Wij(k) in terms of the magnetic 
and electron density power spectra. We define, 

J2«(x) = ((l + 4 1 (y))- 1 (l + (5 n (y + x))- 1 ), 
i? B (x) = (<SB(y)-5B(y + x)>. 

We define Ps(k) and Pg(k) relative to Rb and R$ by the 
relation, 



i?(x) = ^- J d 3 kP(k)t 



2-7rzk-x 



We also define, 

i? N (x) = (Sn c (y)Sn c (y + x)) , 

and corresponding Fourier transform, P^. We separate 
the fluctuations in n e into short and long wavelength fluc- 
tuations, S s and Si respectively, and assume they are ap- 
proximately independent. Thus n c = n e (l + Si) (I + S s ). 
This parameterisation is consistent with the assumption 
that n e /n e has an approximately log-normal distribution. 
We assume that the short wavelength fluctuations are lin- 
ear up to some cut-off scale k^ n . Above this spatial scale 
we have the long wavelength fluctuations, which arc not 
necessarily linear. We then have that over small spatial 
scales, <C k^, 



((l + 5 s (y))- 1 (l + <5 s (y + x))- 1 ), 



\ ll c I 

*(i)(i)V»,- 2 *N(x)], 

= /S[l + f>- 2 B B (x)]. 
It follows that for k kn n we have approximately, 

P s (k)^I 2 n- 2 P N (k). 
We assume that fcii n -C k\, fc 2 . Finally we define, 

PBs(k) =Pbs(V = ^ J d 3 pP S (p)P B (\\k-p\\ 

We note that, 

£ J d 3 kP BS (k)e 2 ^ = (<5B(y) • <5B(y + x)) 
•((l + <5 rl (y))- 1 (l + ( 5„(y + x))- 1 ). 

From the definition of P*°* and Wij(k) we have, 



Wii(k) 



B rcs B rcs 

* 3 



l f 

+ 3% 



qPs(q)dq 
qP B s(q)dq. 
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All that remains is to evaluate 
r 

qP BS (q)dq = / d 3 p 



where, 



3 Ps(p) 



d K 



Pb(v) 



4?r Jv k (p) 47r ii r +pir 

where Vfc(p) is defined by ||r + p|| > k. We approximate 
this integral using ||r+p|| »j)inO<r<p and ||r+p|| w 
r in < p < r, and find 

qP B s(q)dq 



Qo 



2B rcg uj 



V/, 



COS 



(2A) 



Ik 



(SB 2 ) 

+1% 



qP S (q)dq 



POO 

/ qP B (q)dq. 
Jk 



Defining, 

VFb/n 
we have for An, 



= r 
jk 



qdqP B/N (q), 



4 \M eS m 2 eS J 

1 / 2B re gUJ 
4 \M cS fh 2 cS 

We can also calculate the variance of the Qi. We define, 

of (Al A 2 ) - (Q i (A 1 )Q i (A 2 )) - Qi(Ai)Qj(A 2 ), 
and find that, 

4a 2 = \\A XX {A U A 2 )\\ 2 + \\A XX {A U -A 2 )\\ 2 
+ ^(Ar^f + IK^-A,)!! 2 
- ||A y (A 1 ,A 2 )|| 2 -||A,(A 1 ,-A 2 )|| 2 



e 2iA ~ Aij(Ai, A 2 ) 



4j| r B 2 -^iu 2 
2M c 2 ff m4 ff 
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- p ya; (Ai,A 2 



|.4„ x (Ai,-A 2 ) 



4i?r s sr s ^i^ 

■ cos 2 A+ 
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cff TO cff 

B?,j£Lsinc2A_ 

+ - — [Wjv(fci) 



+ 



16M>2 

(5y/^Lsinc2A_ 
48Af 2 ff 



Wjv(fc 2 )] 



[W^(fci) + W B (fc 2 )], 



where B 2 - = B^Bf 



(5y (<5B 2 ) /3. It is now straight- 
forward to calculate the expectations of P 7 „^ and the 
Qi at a single frequency corresponding to A, since 



4a 2 = 2Re(^(A 1 ,A 2 )A /s/ (A 1 ,A 2 )) 

+ 2Re(^(A 1 ,-A 2 M w (A 1 ,-A 2 )) 

+ 2Re(^(A 1 ,A 2 )^(A 1 ,A 2 )) 

+ 2Re(^(A 1 ,-A 2 )^(A 1 ,-A 2 )), 

4a 2 - 2Re(^(A 1 ,A 2 )^ w (A 1 ,A 2 )) 

- 2Re(^(A 1 ,-A 2 )^(A 1 ,-A 2 )) 

- 2Rc(^(A 1 ,A 2 )^(A 1 ,A 2 )) 

+ 2Re(^(A 1 ,-A 2 )^ !/x (A 1 ,-A 2 )). 



2q 



(A XX (A,A)+Ayy{A,A)) 

(A XX (A,A) -Ayy(A,A)) 

(A xy {A,A) + Ay X (A,A)) 
(A xy {A,A)-A yx (A,A)) 



Note that we have assumed that the Ai have an approx- 
imately Gaussian distribution so that the expectation of 
four Ai can be written in terms of the expectation of 
two Ai. If we take Ai = A 2 = A then when P 7 „^, is 
dominated by the terms proportional to Wn, we have 



U2 , 1^2 



-^ + -Q 2 cos40 reg , 



-V 2 , 



^Qo cos40 rcg , 



1 - 



Note that this implies Q v = 0. It follows that, 



2 ' 7-0 



2 - -O 2 



1 ( 2B cS lo 



2 VM cff m 2 ff 
_1 / 2B rcg u \ 
4 ^Meffm 2 ff y 



r 3 



COS 



(2A) 



8m; 



cff«S 



24M 2 ff 



/N^B(fccrit), 



where B c 2 ff ee B r 2 cg /2 + (<5B 2 ) /3, and fc crit = \A\/wL. 
Similarly, 



2q 



Qo cos26> rog , 
Qo sin26> rcg , 



We note that in the cell model we had of w Vj^^/2. 



3. Measurement Issues 

The above evaluation for the cell and power spectrum 
models is appropriate if measurements of the Stokes pa- 
rameters and hence Qi probe only a single frequency at a 
time, i.e. they are precisely monochromatic. In practice, 
a measurement at a nominal frequency u>, will generally 
be an average over some frequency band with width 5w. 
We call 8u)/u> the spectral resolution. We define 2<5A to 
be the difference in the value of A across the frequency 



24 



band. When measurements of the Stokes parameters are 
made over a finite frequency band, the observable vari- 
ance in the Qi is no longer given by the of expressions 
presented above but, in the (more general) case of the 
power spectrum model, by 

i pSA pSA 

(ZOlA) J_ SA J-SA 



We find that for both the cell model and power spectrum 
model the largest term in of when <5A = is proportional 
to l/d A when SA > 1. If SA > 1, then a\ < V*^. 
Similarly the variance of 'P 7< -^0 is also much smaller than 
~P%<^tb when SA > 1 so that, to a good approximation, 
we have V-y^c/, — P 7< _0. 
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